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Abstract We evaluate the current status of supernova remnants as the sources 
of Galactic cosmic rays. We summarize observations of supernova remnants, 
covering the whole electromagnetic spectrum and describe what these obser- 
vations tell us about the acceleration processes by high Mach number shock 
fronts. We discuss the shock modification by cosmic rays, the shape and 
maximum energy of the cosmic-ray spectrum and the total energy budget of 
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cosmic rays in and surrounding supernova remnants. Additionally, we discuss 
problems with supernova remnants as main sources of Galactic cosmic rays, 
as well as alternative sources. 

Keywords Supernova remnants • cosmic rays • acceleration of particles 



1 Introduction 

This year, we celebrate the discovery of cosmic rays (CRs) a century ago. In 
1912, Victor Hess went up in the atmosphere with a balloon, and found that 



the strength of ionizing radiation increased as he went up (Hess 1912). He 
realized that this increase points towards a cosmic origin of the radiation. 
For this reason the radiation was dubbed 'cosmic rays'. Later, it was found 
that these CRs were in fact energetic particles, rather than electromagnetic 
radiation. However, 'cosmic rays' is, although not an adequate description, 
still the currently used name for these particles. 

Since their discovery, CRs have been subject to intensive study. For a 
good reason, as CRs can have tremendous energies: particles with energies 
more than 10 20 eV have been detected (Figure [l]). Additionally, their energy 
density in the interstellar medium is about one third of the total energy 
density in the interstellar medium, comparable to the thermal and magnetic 
energy densities. Moreover, CRs provide an additional source of observational 
information on their accelerators next to electromagnetic radiation. 

Supernova remnants have been the main candidates for accelerating Galac- 



tic CRs for decades (Baade and Zwicky 1934), but this is still subject to 



debate. In this paper, we describe the current observational status of CRs 
in SNRs. We describe the emission directly originating from the CRs as well 
as how we can get information on CR acceleration in remnants from other 
observational characteristics. Also, we briefly discuss the problems with su- 
pernova remnants as main sources of Galactic CRs as well as alternative 
sources. 



1.1 Supernovae and their remnants as sources of Galactic cosmic rays 

Particles with energies less than 10 18 eV have gyro radii smaller than the 
thickness of the Galactic disk. Hence, particles with higher energies are likely 
to escape the Galaxy, indicating that the CRs well above 10 18 eV are prob- 
ably accelerated outside of the Galaxy. In contrast, the particles with lower 
energies are contained within the Galactic magnetic field. The trajectories of 
these particles are curved and randomized by the magnetic field, causing their 
direction of arrival at Earth to be largely uncorrelated with the direction to 
their sources. 

One clue towards the sources of Galactic CRs is the shape of the CR 
spectrum. The constant power-law slope of 2.7 over many decades in energy 
likely indicates that the particles of different energies were accelerated by a 
similar mechanism. As this is a steep power law, most of the energy in CRs is 
contributed by CRs with low (GeV) energies. Unfortunately, the spectrum at 
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these lower energies is heavily modified by solar modulation. Therefore, it is 
hard to determine the total energy budget in CRs from experiments directly. 
Nevertheless, | Webber ( 1998 ) reports an CR energy density in the Galaxy of 
1.8 eV per cubic centimeter, correcting for solar modulation effects by using 
Voyager and Pioneer data. Isotope studies show that CRs with GeV energies 
diffuse out of the Galaxy in 15. 0± 1.6 Myr (Yanasak et al. 20011. Galactic 



supernova explosions can replenish these losses, provided that they transfe r 
- 10-20% of their kinetic energy (£ S n = 10 51 erg) into CRs (Hillas] |2005|>, 



assuming two to three supernova explosions per century (Tammann et al. 
19941). 



During a supernova explosion, the ejecta are expelled into the surrounding 
medium with speeds as large as tens of thousands of kilometers per second. 
The ejecta carry the kinetic energy of the supernova explosion and cause 
high Mach number shock fronts in the interstellar medium. At the end of the 
seventies, four independent papers proposed an acceleration mechanism for 
charged particles that operates at high Mach number shocks ( Axford et al 



1977 Krymskii 1977 Bell 1978 Blandford and Ostriker 1978 ). This mech- 



anism concerns particles scattering up and down the shock front on turbulent 
magnetic fields, thereby effectively gain energy each time they cross the shock 
front. This process, called 'diffusive shock acceleration', natura lly produces 
energetic particles, with a power-law spectrum with index 2 (see Malkov and 
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Fig. 1 Cosmic-ray spectrum as observed by different experiments. Data compiled 
by J. Swordy, courtesy to Klara Schure. 
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Table 1 Table with shell-type supernova remnants that hav e X-ray synchrotron 
e mission, together with t heir properties a t othe r wavelenths. 1 Abdo et al.|(|2010b|), 
2 |Albert et al.| J2007b||; |Acciari et al.j J2010|), ^Fesen et al. (|2006|), 4 |Reynolds| 
let al.| fl2008|>, a |H.E.S.S. Collaboration et all fl201lD, |Aharonian et al.| (2009| 7 



Stephenson and Gree n (2002), 



Abdo and Fermi LAT Collaboration 



(2011) 



Aharonian et al. (2004|), lu |Wang et al. (1997|), 11 Tanaka et al. 


(201l|, ^|hess 


Collaboration: F. A. Aharonian| (|2006|), " Katsuda et al.| (|2008 


i, 14 Acero et al. 



Drury 2001 Schure et al. 2012 for a review). The somewhat softer CR 
spectrum observed in the Earth's, might be explained by propagation effects 
of CRs through the Galaxy. 

As the direction of arrival of Galactic CR is essentially uncorrelated to 
their sources, most of the observations on the sources of Galactic CRs, comes 
from electromagnetic radiation. As energetic electrons gyrate in the local 
magnetic field, they emit synchrotron emission. This synchrotron emission 
is observed at radio wavelengths in supernova remnants, revealing the pres- 
ence of GeV electrons ( Sh klovskiij |1953| |Minkowski| 1957). Also, X-ray syn- 
chrotron emitting TeV electrons were identified at the shock fronts of the SN 
1006 supernova remnant and others (Koyama et al. 1995 and see section 



2.3 ). In the last decade, GeV and TeV 7-ray emission has been observed from 



several supernova remnants (Table |1.1|) , indicating 


; the presence of particles 


with GeV and TeV energies (e.g., Aharonian et al. 


2004| 


Abdo et al. 


2010b 



and see section 2.13). 



1.2 Different types of supernova explosions 



Classification of supernova explosions is traditionally based on the presence 
of hydrogen lines in the optical spectrum at maximum light. A supernova 
with hydrogen lines in its spectrum is classified as a type II, and without is 
classified as a type I (Minkowski 1941 ). More specifically, a spectrum without 



hydrogen lines but with strong Si II lines is classified as a type la spectrum, 
without Si II lines but with strong helium I lines as type lb and if the spectra 
are lacking both strong Si II and He I lines, the supernova is classified as a 
type Ic. Type II supernovae have subclassifications as well, depending on 
their light curves. A supernova explosion with a plateau in its light curve is a 
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type IIP, one with a linear phase is a type IIL (Barbon et al. 1979). Another 



subcategory is type lib supernovae: their spectrum shows a weak hydrogen 
line shortly after the explosion, which is absent in the second 



maximum 



of the lightcurve (Filippenko 1988). Recently, spectra of light echoes of the 



explosion that resulted in the Cassiopeia A remnant (Cas A) showed that this 
was a type lib supernova ( Krause et al. 2008 ) . A fourth category of supernova 



explosions are type Iln supernovae ( Schlegel 



1990). Type Iln supernovae are 



characterized by narrow emission lines in their spectra, probably produced 
by the interaction of the ejecta with dense circumstellar medium around 
these explosions, expelled by the progenitor star before explosion. For a more 
extensive description of observational characteristics of supernova explosions, 
Filippenko| (|1997|). About 24% of all supernovae are type la's, 10% 



sec 



<IYV 



type Ibc's and 57% are type IPs (Li et al. 2011) 



Apart from this observational classification scheme, we now know that 
type la supernovae are related to thermonuclear explosions of white dwarfs 
(first suggested by Hoyle and Fowler 1960). White dwarfs explode when 



they reach the Chandrasekhar mass of 1.38 Mq by accreting material from a 
binary companion. Although the link from type la supernovae to thermonu- 
clear explosions of white dwarfs is fully canonized, it is rather unclear which 
progenitor systems (single or double degenerate binaries) create these explo- 
sions. More specifically, it is unclear which systems accrete efficient enough 



to explain the o bserved frequency of type la explosions (see Hillebrandt and 
Niemeyer |2000 for a review). 



The other (non type la) explosions mark the deaths of massive (M > 
8M Q ) stars. These are called core collapse supernovae. After these stars have 
burned all the nuclear energy in their cores, they are left with an iron core. 
The gravitational pressure will result in an implosion of this iron core to a 
neutron star or a black hole. The envelope will bounce onto the core into the 
circumstellar medium with tens of thousands kilometers per second. 

Although the explosion mechanisms for core collapse and thermonuclear 
explosions are quite different, the kinetic energies dumped into the ambient 
medium are reasonably similar: ~ 10 51 erg. 



1.3 Supernova remnant evolution 

After a star goes supernova, the stellar ejecta expand into the ambient 
medium. These expanding ejecta will create a shock front, the forward shock, 
which will slow down, as it sweeps up an increasing amount of ambient 
medium. As a reaction, a second shock will develop, moving away from the 
forward shock, into the stellar ejecta. This second shock is called the reverse 
shock. The evolution of supernova remnants is characterized by the free ex- 
pansion phase, the Sedov phase and the radiative phase (in chronological 
order). Initially, the forward shock velocity is approximately constant as long 
as the swept-up mass is smaller than the ejecta mass (free expansion phase). 
Precisely speaking, the forward shock velocity slightly decelerates during the 
free expansion phase because it transfers kinetic energy into sweeping up the 
ambient medium (IChevalier 1982b). In the subsequent phase (Sedov phase) 



most energy has transferred from free expanding ejecta to the shock-heated 
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shell. As radiative cooling is still negligible in this phase, the supernova rem- 
nant expands adiabatically in this phase. As the shock slows down to veloci- 
ties below 200 km s^ 1 , the post-shock temperature will be less than 5 x 10 5 
K. This will cause the radiation of H, He, C, N and O line emission to increase 



drastically ( |Raymond| |1979| |Draine and McKee[ |1993| |Schure et al| 2009[ ), 
and therewith radiative cooling becomes significant (radiative phase). The 
length of each evolution phase of supernova remnants depends on supernova 
type and structure of their surroundings. 

For type la supernovae, the surrounding medium typically has a constant 
density ISM with no ~ 1 cm -3 and To ~ 1 eV, where no and To are the 
number density and the temperature of ISM, respectively. A mass M j ~ 
1 Mq ejected with E$n ~ 10 51 erg in a supernova explosion has an initial 
velocity (V ) of 



Pi \Uf* erg) \1M B J 

The forward-shock velocity during the free expansion phase is V s ~ Vq, so 
that the Mach number is given by 

V, ,„ 3 ( -Esn \"7 M, s \""V To ^- 1 ' 2 



where c s = 10 6 cm s _1 (To/1 eV) is the sound speed. The ejecta start to 
decelerate when the swept-up mass becomes comparable to the 

ejecta mass, where R s and po are the forward shock radius and the density 
of the ambient medium, respectively. The transition from the free expansion 
phase to the Sedov phase occurs at shock radius -Rscdov, 

and at time fsedov, 

flsedov 01 „ ( E SN y 1/2 ( M cj V /6 / n x-i/3 

(4) 



For supernova remnants in superbubbles, i?s e dov and ^Sedov are about ten 
times larger than that in a typical ISM because the density is no ~ 10~ 3 cm~ 3 . 

After tscdov, the supernova remnant evolution is in the Sedov phase as 
long as the radiative cooling is negligible. Assuming that the swept-up mass 
is larger than the ejecta mass and the total explosion energy is conserved, 

\ v * x y p ° r3s = const ' ' (5) 

the evolution of the forward shock, R s and V s are given by 

( 1 v /5 

Rs = ^Sodov I ~ , (6) 

V iSodov / 
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v s 



dR s 
dt 



2i? s , 



cdov 



5ts. 



cdov 



^Scdo 



-3/5 



(7) 



A simple treatment of the dynamics of the for ward shock from the f ree ex - 
pansion to the Sedov phase has be en gi ven by Ostriker and McKee (1988), 
Bisnovatyi-Kogan and Silichl (|1995j) andlTruelove and McKcc (1999T. 



When the supernova remnant age becomes comparable to the radia- 
tive cooling time, the supernova remnant evolution changes from the Se- 
dov phase to the radiative phase (the pressure driven snowplow phase and 
the momentum-conserving snowplow phase). This transition occurs at t tI , is 
given by (|Petruk[ |2005| 



2.8 x 10 4 yr 



/ n \ 4/17 

V10 51 ergy VI cm- 3 / 



-9/17 



(8) 



The asympto tic supernova remnan t evolu tion of the radiative phase was in- 
vestigated by McKee and Ostriker (1977) (R s oc t 2 ^ 7 for the pressure driven 



snowplow phase) and Oort ( 1951 ) (R s oc i 1 / 4 for the momentum-conserving 



snowplow phase). [Bandiera and Petruk ( 2004 1 found more accurate analytical 
solutions for the radiative phase. Truelove and McKee ( 1999 ) also provided 



approximate solutions for the reverse shock evolution. For uniform ejecta and 
a uniform ISM, the reverse shock position, i? r , and the velocity of t he reverse 
shock in the rest frame of the unshocked ejecta, V^j is given by (Truelove 
and McKeel 119991) 



Ri — Rs~ 



X 



1 tScd, 

1.29 <! 1 + 0.947 



0.707 1.22 - 0.0465 



0.533 In 



*Scdo 



(t < iSedov) 
(t > iSedov) 



(9) 



I!., 



X 



*Sod. 

1.22 
0.377 



L 



*Sod 

-1-0.0329 



) 3/2 {l + 0.947 (-^)" 



■5/3 



Sedov / 



(t < isedov) 
(t > tsedov) 



(10) 



The velocity of the reverse shock in the rest frame of the unshocked ejecta, V x , 
increases with time for uniform ejecta. Truelove and McKee ( 1999 ) provided 
approxim ate solutions for some cases of type la supernovae evolution (see 
table 7 oflTruelove and McKeel 119991). 



1 Note that this is the velocity relevant for the physical conditions behind the 
shock front (outside of the reverse shock) . 
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Fig. 2 Evolution of the maximum energy. The black, blue and green lines show 
the age-limited, the cooling-limited and the escape-limited maximum energy, re- 
spectively. The red dashed lines show the evolution of the maximum energy of 
CR electrons. Thick and thin lines show cases where the cooling-limited maximum 
energy increases or decreases with time during the Sedov phase, respectively. 



For core-collapse supernovae, the structure of surrounding medium is 
more complicated and depends on the details of the stellar mass loss his- 
tory, which in itself depends on the initial mass, binarity and local ambient 
medium density (e.g. van Marie et al. 2004 1. Roughly speaking, in order of 



increasing distance to the center, the ambient material is organized by the 
dense red supergiant wind, a shocked low density wind created during the 
main sequence stage, and interstellar medium. The density structure in the 
dense red supergiant wind is po oc r~ 2 whereas density structures of other 
regions are constant. Approximate solutions of supernova remnant evolution 
in these media are also obtained in the same manner as type la supernovae 



(e.g. 


Chevalier 


1982a 


Truelove and McKee 1999 


Ostriker and McKee 


1988 


Bisnovatyi-Kogan and Silich 


1995 


Laming and Hwang 


2003 


1. 



1.4 Cosmic-ray acceleration efficiency during supernova remnant evolution 

As mentioned above, there are three phases for supernova remnant evolution. 
Although we do not understand when particle acceleration stops, observa- 
tions tell us that relativistic electrons are accelerated during the free expan- 
sion phase and the Sedov phase. Particle acceleration mechanisms during the 
radiative phase have not been investigated in detail. However, these are not 
expected to provide a significant contribution to the Galactic CR density, as 
the shock velocity is low, and a substantial fraction of the shock energy is 
lost by radiative losses. 

The maximum energy of CR protons is limited by the finite age or the 
escape condition. According to the theory of diffusive shock acceleration, 
an estimate for the typical time to accelerate a particle with an elementary 



1983 


Jokipii 


1987 


Parizot 
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etaL||2006| : 

■ '"^ ^ =124 Vg Bzl / V ' E 



yr, 



' V s 2 r-l ls -HsOOOkms- 1 / V100TeV/r 4 -i 

(11) 

with D2 the downstream diffusion coefficient for particles with energy E 



(Drury 19831, V s the shock velocity, B_4 the downstream magnetic field in 
units of 100 /iG, r) g is the gyro factor and r 4 the overall compression ratio 
in units of 4. The factor 1.8 comes from taking into account the difference 
between magnetic field upstream and downstream of the shock. Here, we as- 
sume that there is a highly turbulent, isotropic, magnetic field upstream and 
that the downstream magnetic field strength is determined by compression 
of the perpendicular component of the upstream magnetic field. Note that a 
higher compression ratio corresponds to a longer acceleration time, but that 
a high compression ratio itself is the result of efficient acceleration. 

To get some understanding on the typical properties of particle accelera- 
tion we can approximate equation |11| by: 

_ D _ cE 

^acc Vacc ^ 2 ^ aC °^ S 3ei3V^? ' 

where ?7 acc is a numerical factor which depends on the shock compression 



ratio and D is the diffusion coefficient around the shock (Drury 1983). The 
diffusion coefficient relates to the mean free path (A) as follows: D = ^Ac, 
where c is the speed of light. We generally express A in terms of gyroradius 
(r g ): A = TgVg- The lower r/ gl the more turbulent the magnetic field. r/ g = 1 
is generally referred to as Bohm diffusion. 

One can obtain the age-limited maximum energy, -Em^ge, from the con- 
dition t acc = t assuming r& = 1 

_ 3eBV s H 

77acc%C 

The escape time due to the diffusion, t csc , is written by 
D rj g cE 

where ry osc is a numerical factor (ry csc < 1). Note that t esc is similar to t age 
for ^?7csc?7acc ~ 1 an d R s /V s t s» 1. 

For this reason, the escape-limited maximum energy is similar to the 



age-limited maximum energy (Drury 2011). The time dependence of the 



maximum energy B mjage « S m>esc is given by 
_ B(t)V B (t)H 

where we assume that ry acc is constant. For a constant density profile, the am- 
plified magnetic field and their fluctuations are thought to be constant during 
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the free expansion phase because the shock velocity is approximately con- 
stant. Therefore, the maximum energy linearly increases with time, E m ^ ge oc 
t. After the free expansion phase (t > iscdov), the shock velocity decreases 
with time, so that the maximum energy decreases with time (E m ^ sc oc 
BVgt oc Bt~ x l b ). Additionally, the magne tic field is thought to decrease 
with time (B oc or oc V^ 3 , see section 2.7), leading to an even faster 



decrease of the maximum energy. 

The number of CRs accelerated by supernova remnants increases with 
time since the swept-up mass increases with time. However, the detailed 
evolution of the number of CRs is unclear. Recently, several authors discuss 
the spectrum of escaping CRs corresponding to the source of Galactic CRs 



(Ptuskin and Zirakashvili 2005 Ohira et al. 2010 Caprioli et al. 2010 



Drury 



201 ip . They showed that the spectrum of escaping CRs depends on 



the evo. 

moreover, (iPtuskin and Zirakashvili 



ution of the number of CRs and the evolution of the maximum energy, 



2005 ) showed that the maximum energy 



to which particles can be accelerated differ for different types of supernovae. 



2 Observational diagnostics 

2.1 Shock parameters 

As mentioned before, in order to maintain a steady Galactic CR energy 
density, on average 10-20% of the explosion energy of a supernova has to 
be converted to CR energy. For these efficiencies, the acceleration process 
changes the shock structure and makes a precursor in the upstream regiorF] 



(Drury and Voelk[|1981|). In the precursor region, magnetic-field amplification 


(Lucek and Bell[ |2000 Bell, 2004 Ohira and Takahara 




2010 


Bykov et al.[ 


2011b Schure and Bell 201lTReville and Bell[ |2011|) 


and plasma heating 


(McKenzie and V61k[ |1982| Ghavamian et al. 2007; Rakowski et al.| 2008 


Niemiec et al.j 2008 |Ohira et al. 


2009a| Raymond et al. 


2011) are thought 



to occur. Although there is more and more evidence that magnetic field 
amplification is important, the process itself is not yet fully understood. 



In the shock rest frame, the upstream plasma is decelerated by the CR 
pressure and and heated by plasma waves at the precursor region. As a 
result, the effective Mach number becomes small at the jump of the shock 
(subshock) and the compression ratio, at the subshock becomes less than 
four. On the other hand, the total compression ratio (the density ratio of the 
far upstream to the downstream) becomes large (r > 4). This is because a 
part of the downstream pressure is due a) to the relativistic particle with the 
adiabatic index of 7 = 4/3 (r = 1 + 1/7 - 1) and b) escaping CRs take away 
shock energy. 

The energy spectrum of accelerated particles also changes from that pre- 
dicted by standard diffusive shock acceleration. The energy spectrum of low 
energy CRs becomes steeper than that of the standard diffusive shock accel- 
eration because low energy CRs sample only the contrast in velocity across 

2 The upstream region is the region in front of the shock, i.e. outside the remnant 
for the forward shock. 
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the subshock (r < 4). On the other hand, the spectrum of high energy CRs 
becomes harder than that of the standard diffusive shock acceleration be- 
cause high energy CRs sample the total velocity difference in total shock 
region (subshock plus shock precur sor: r > 4) because of the longer diffusion 
length (Berezhko and Ellison 1999 ). As a result, the energy spectrum in wide 



energy range is not a single power 



Malkov and Drury 2001 and section 2.11 ) 



aw form (e.g. Berezhko and Ellison 1999 



If a large fraction of the shock energy is converted to the CR energy, the 
mean temperature of the downstream plasma becomes smaller than that of 
the standard shock, (T) = 3{j,m p Vs /16, where \i is the mean particle mass 
(see section 3.3 ). 



2.2 Radio observations 



So far, about 300 supernova remnants have been observed in the radio band 
( Green 2009 ) . Their spectra are nonthermal and their radio spectral indices^ 
a [S v <x v~ a ) 7 are typically a = 0.5 ± 0.1 (see figure 1 of | Reynolds et aF 



2011). Radio spectra can be interpreted as synchrotron radiation from rel- 
ativist ic electrons accelerated by diffusive shock acceleration. |van der Laan| 
( 1962 ) considered the adiabatic acceleration of ambient CR electrons by the 



radiative shock and the compressed magnetic field. The spectral index pre- 
dicted by the van der Laan mechanism depends on the spectral index of low 
energy CR electrons with energies around 30 MeV (Cox et al. 19991. How- 



ever, spectra of low energy ambient CR electrons have not been understood 
completely. The radio luminosity and spectra of young supernova remnants 
cannot be understood by the van der Laan mechanism, and requires active 
acceleration of electrons by the supernova remnant shocks. 

Time evolution of radio fluxes is observed in some young supernova rem- 
nants. The radi o flux of Cas A has been obse rved for 50 yr and decreases 
with time (e.g. Helmboldt and Kassim 20091. This means that either the 



magnetic field or the number of CR electrons decreases with time. Not only 
the evolution of radio flux but also the evolution of the radio size is observed 



for very young supernova remnants (e.g. Weiler et al. 2002 Bartel 2009 



Ng et al. 2008). Therefore, we can compare evolution models of supernova 



2011) 



remnants and observations (e.g. Marti- Vidal et al.| 

Polarization of radio synchrotron radiation has been also observed. Young 
supernova remnants have a predominantly radial magnetic field structure, 
whereas old supernova remnants have a circumferential magnetic field struc- 
ture (Dickel and Milne 1976 Milne 1987). The circumferential magnetic 



field structure of old supernova remnants can be interpreted as th e com- 
pression of the upstream magnetic field by the radiative shock. Gull ( 1973 ) 



suggested that the radial magnetic field of you ng supernova remnants re sults 
from stretch by the Rayleigh- Taylor instability. Jun and Norman ( 1996 ) con- 



firmed that the Rayleigh- Taylor instability radially stretches the magnetic 



3 Note that the power-law index of a radio spectrum (a) is defined for a flux 
energy density, whereas in X-rays one use the spectral index (_T) for the number 
density flux, the relation between the two indices is _T = a + 1. If the number 
density of the electron distribution is labeled q, we have that q — 2a + 1 = 2_T — 1. 
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field. However, Schure et al. (2010) showed that the Rayleigh- Taylor insta- 
bility in the gas with the adiabatic index of 5/3 can not explain the radial 
field polarization observed near the forward shock, whereas the case of the 
adiabatic index of 1.1 can reproduce that. In order to make the adiabatic in- 
dex of 1.1 in young supernova remnants, escaping CRs have to carry away a 
substantial fraction of the shock energy. |Zirakashvili and Ptusk in ( 2008 ) pro- 
posed an alternative scenario that upstream density fluctuations produced 
by the non-resonant instability radially stretch the amplified magnetic field 
in the downstream region. 



2.3 X-ray observations 



Early observational evidence that supernova remnants are capable of accel- 
erating particles to high energies mainly consisted of the radio synchrotron 
radiation from relativistic electrons from most supernova remnants. But the 
electrons responsible for the radio emission have energies in the order of a 
GeV. The discovery by the ASCA satellite that the featureless X-ray spec- 
trum of SN 1006 is caused by synchrotron radiation from 10-100 TeV elec- 



trons (Koyama et al. 1995) suddenly increased the range of determined elec- 



tron energies by orders of magnitude. 

The X-ray emission from SN 1006 was long recognized to be peculiar, as 
its overall spectrum appeared to be devoid of line emission. Indeed , a syn- 
chrotron interpretation was offered by Reynolds and Chevalier ( 1981 ). But a 



sophisticated th ermal explanation, invo lving fully ionized carbon-rich super- 
nova ejecta, by Hamilton et al. ( |1986 ) was regarded as the most likely ex- 
planation. ASCA, however, revealed that the X-ray emission was dominated 
by continuum emission from near the shock front. Hot supernova ejecta are 
present in the interior of the supernova remnant, and produce, among others, 
oxygen and silicon line emission, which was hard to detect in the synchrotron 
dominated spectrum from the supernova remnant as a whole. The presence 
of oxygen-rich ejecta was clearly in contradiction to the model of Hamilton 
et al. (1986). Moreover, the ASCA observations showed that the the contin- 



uum emission was associated with regions close to the forward shock front; 
an unmistakable signature of synchrotron radiation from electrons that are 
being actively accelerated by the supernova remnant shock. 

Around the same time as the discovery of X-ray synchrotron emission 
from SN 1006, it became clear that some young supernova remnants emit 
non-thermal, hard X-ray emission. This was in particular apparent for Cas 
A, which is detected in hard X-rays from ~ 10 — 100 keV, indicating a nearly 
power-law like spectrum with index ~ 3.2 ( The et ah] |1996[ |Allen et al. 



1997 Favata et al. 1997 Vink et al. 2001 Renaud et a 



2006). For other 



supe rnova remnants h ard X-ray emission up to about 30 keV was detected 
(e.g. Allen et al. 1999 ). Although a non-thermal bremsstrahlung explanation 
for the hard X-ray emission of some supernova remnants exists (e.g. Asvarov 



et al. 1990 Favata et al 



1997 Laming 2001 Bleeker et al. 20011, the 



rapid energy exchange of the ~ 10 — 100 keV electrons that cause X-ray 
bremsstrahlung with background electrons, makes that a power-law shaped 
bremsstrahlung spectrum is unlikely (Vink 2008b). 
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Since 1995, the number of X-ray synchrotron emitting supernova rem- 
nants has grown rapidly. Among the earliest identifications of X-ray syn- 



chrotron emitting supernova remnants were RX J1713.7-3946 (Koyama et al. 



1997 Slane et al. 1999) and RX J0852.0-4622 (Slane et al. 



200ip , two su- 



pernova remnants that later turned out to be bright in TeV gamma-rays 
(Aharonian et al. 2004 2005). Also the large diameter supernova r e mnan t 
RCW 86 showed signs of non-thermal X-ray emission flVink et al.j [1997]) 



which was later argued to be synchrotron radiation (Bamba et al. 2000 



Borkowski et al. 2001) 



The angular resolution of Chandra was need ed to find that also the 
youngest supe rnova remnants l ike Ty cho/SN 1572 (Hwan g et al. 2001 1, Ke- 
pler/SN 1604 ([Reynolds et al.||2007) and Cas A ( |Gotthelf et a^|2lM[|Vmk 



and Laming 2003 ) emit X-ray synchrotron emission. In those supernova rem- 



nants this non-thermal emission is concentrated in very narrow filaments 
(1-3") marking their shock fronts. 



2.4 Conditions for X-ray synchrotron emission 

Apart from proving that supernova remnants can accelerate particles up to at 
least 10 14 eV, X-ray synchrotron emission also provides a powerful diagnostic 
for measuring the acceleration conditions. The presence of X-ray synchrotron 
filaments near shock fronts suggests that the electrons have been accelerated 
during the last century, or even during the last decade. This identifies su- 
pernova remnants as objects that are actively accelerating particles, rather 
than just objects that contain particles accelerated during an earlier phase. 
Moreover, as we shall see below, X-ray synchrotron emission is only possible 
if acceleration is surprisingly fast, with a diffusion coefficient that is rather 
small, and within an order of magnitude of the Bohm limit. To understand 
the conditions under which X-ray synchrotron emission can occur, we need a 
closer look at the acceleration properties in supernova remnant shocks, and 
at the characteristics of synchrotron emission. 



Equation 11 describes the time it takes (t aC c) to accelerate a particle 
to an energy E. The maximum energy an electron (or proton) can obtain 
is limited by the age of the shock. For example, for a supernova remnant 
with _B_4 = 0.1, V s = 5000 km s _1 an electron (or proton) can only be 
accelerated in ~ 1000 yr to 100 TeV. If the maximum electron energy, and 
the corresponding synchrotron cut-off frequency, is indeed determined only 
by the time available for acceleration, the synchrotron spectrum is said to be 



limited (Reynolds 1998) 



However, unlike protons, electrons lose their energies quite rapidly through 
synchrotron radiation and inverse Compton scattering, certainly for energies 
in excess of 1 TeV. So another limit on the maximum energy of electrons is 
set by the condition t acc rj t\ oss . If the maximum energy for electrons is set by 
this condition the synchrotron cut-off frequency is said to be loss-limited. In 
that case the ions can, in principle, obtain much higher energies. For magnetic 
fields in excess of ~ 3/iG synchrotron losses usually dominate over inverse 
Compton losses, unless the local radiation field is in strong excess over the 
microwave background radiation field. 
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The synchrotron loss time for an electron is given by: 
E ,/ E 



12.5 BZ\ 



( * V 

V100 TeV/ 



yr. 



(16) 



y dE/dt " ~*Vl00TeV, 
The cut-off photon energy for the loss limited case can be calculated by 



combining Eq. 11 and Eq. 16 which gives 



w ™ TeV - o.3 2 , ! -( i ^ s )-(^ r ) (^iy<\ 



V 100 fiG 



(17) 



where is the compression ratio divided by 4. This equation can be combined 
with the relation between the peak in synchrotron emissivity, the so-called 
characteristic frequency f c h, for an electron with given energy E moving in 
a magnetic field B (e.g. Ginzburg and Syrovatskij 19671: 



v ch = 1.8 x lO 18 ^ 



E 
1 erg 



Hz. 



hp, 



'oh 



13.9 



B± 



E 



100 nGJ V 100 TeV 



keV, 



(18) 



with B±_ ss \j2jZB the magnetic field component perpendicular to the mo- 
tion of the electron. Taken together Eq. [18] and Eq. [17| indicate that the 



synchrotron cut-off frequency is expecte d to be (c.f. Aharonian and Atoyan 



1999 Zirakashvili and Aharonian 2007): 



h v. 



cut— off 



1.471 



V, 



5000 kms- 



kcV. 



(19) 



This shows that in the loss-limited case the cut-off frequency does not depend 
on the magnetic field strength, and that for typical young supernova remnants 
like Cas A, Kepler and Tycho, which have V s ~ 5000 km s -1 , one expects 
to find X-ray synchrotron emission provided that r) s ~ 1. One should realize, 
however, that the the synchrotron emission spectrum for a given electron 
energy is rather broad. As a result a rather sharp, exponential, cut-off of 
the electron spectrum results in a more gradual turn-over of the resulting 
synchrotron spectrum. Nevertheless, all observed X-ray synchrotron spectra 
of supernova remnants are steeper than the radio spectral indices. In X-rays 
the observed spectral indices are r = 2 — 3.5, whereas in radio r = a + 1 = 
1.4 — 1.8. This indicates that in X-rays the spectra are observed near or 
beyond the cut-off frequency. Since for ry g 3> 10 we would not be able to 
detect X-ray synchrotron emission, the detection of X-ray synchrotron has 
already significant consequences: it suggests that the magnetic fields near the 
shocks of young supernova remnants are highly turbulent. 
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0.001 



Fig. 3 Left: Chandra X-ray image of Cas A, showing in green the synchrotron dom- 
inated continuum emission. Note the narrow filaments alineating the shock front. 
Right: The maximum cosmic-ray electron energy versus magnetic field strength for 
the region just do wnstream of Cas A ' s shoc k front, as determined from the thickness 
of the filaments (Vink and Laming 20031. The shaded area is excluded, because 
the filament width cannot b e smaller than the minimum possible diffusion leng th 
(c.f. |Vink and Lami ng 2003, these figures were published before in |Vink[ |2006[ )). 



2.5 The widths of X-ray synchrotron filaments 



The narrowness of the synchrotron filaments, in particular for the supernova 
remnants Cas A, Tycho and Kepler, is very striking. For these supernova 
remnants the width is typically 1-3", which corresponds to typical physical 
widths of ~ 5 x 10 16 cm. It was soon recognized that these widths indicate 



high magnetic fields near the shock front ( 


Vink and Laming 


2003 Bcrczhko 


et al. 


2003a| 


Bamba et al. 


2004 


Ballet 


2006p. One can 


understand this 



using two different approaches. One approach is to consider that the plasma 
containing the highest energy electrons is moving away from the shock front 
with a relative velocity Av — V s /r, with r the shock compression ratio. As 
the electrons are moving outside the reach of the shock front they do not 
gain, but only lose energy through synchrotron radiation/inverse Compton 
scattering. After a certain time they have lost so much energy that they no 
longer emit X-ray synchrotron emission. The width thus corresponds to the 
advection length scale: 



kdv = = 4.9 x 10 lo r 



1 R -2 
4 U -\ 



5000 km/s/ VlOOTeV 



E 



cm, 



(20) 



where we have made use of Eq. 16 In order to put a more certain limit on 
the magnetic field one can eliminate E by using Eq. [18] to obtain: 



adv 



1.8 x 10 17 7 



4^-4 



3/2 



V* 



hv \-l/2 



-)( — ) 



(21) 



,5000 km/s, 

This procedure is graphically illustrated in Fig. [3j A filament width 5 x 
10 16 cm observed at 5 keV, therefore, corresponds to magnetic fields of 100- 



16 



200 fiG, or even more if one realizes that the measured filament width over- 



estimate the physical width due to projection effects (c.f. Berezhko et al 



2003b). This means that for Cas A and Tycho's supernova remnant the mag- 
netic fields are typically 100 — 500 fiG, well above the compressed magnetic 
field of the Galaxy, which would be around 15 /zG. 

The other way to estimate the magnetic field is by noting that, as a result 
of diffusion, a population of particles will typically spread out over a region 
equal to the diffusion length scale. This can be expressed as 



Miff 



i ±_ 

'AV a 



D 2 r 



5.3 x 10 17 Bll?7 g r4 



V s 



5000 kms 



_rV E ) 

V VlOOTeV/ 



(22) 



Here, D 2 is the diffusion coefficient behind the shock front. This puts at 
least a minimum to the possible width of the filaments (see the shaded area 
in Fig. [3|. Interestingly, this methods has a difference dependence on r, E 
and B. Nevertheless, using this estimate gives virtually the same magnetic 
field values as the advection method (see for example |Ballet 2006 Parizot 
eFaLl[2006l). 



In Eq. 22 one can replace E with hv using Eq. 18 (and ignoring the subtle 



distinction between B± and B, 
region anyway): 



which is probably small in the post-shock 



Z diff pa 1.4 x 10 17 S_4 /2 % r 4 



5000 kms" 



hv \V2 
Ikevi Cm ' 



(23) 



This still relies on the idea that the observed photon energies hv cor- 
respond roughly to the characteristic photon energies hv c \ 11 which may not 
be necessarily true given the broad synchrotron emissivity function for a 
given electron energy. However, one can make use of another property of 
loss-limited synchrotron emission. Eq. [11] and Eq. [22] together give 



2.8 



diff 



(24) 



Combining this with the condition t\ c 
Zdiff ~ 1.5(r 4 - -)^adv, 



tare we see that 



(25) 



This assumes that synchrotron losses are dominated by the losses in the 
downstream region, where the magnetic field is higher than in the shock 

4 



2005 Parizot et al. 2006 Vink et al. 2006). For r 



precursor (c.f. Vink 
this gives /diff ~ (adv This shows that estimating magnetic fields through 
the advection length scale is almost equivalent to estimating them using the 
diffusion length scale. 

Using Eq.[20]and[22]to eliminate E one finds an expression for calculating 
the typical downstream magnetic field that is independent of V s and only 
weakly dependent on the compression ratio and rj g . But it does rely on the 
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assumption that the electrons causing the radiation have energies around the 



cut-off energy (e.g. Parizot et al. 2006 Vink et al. 2006) 



B 9 



26 



adv 



1.0 x 10 18 cm/ 



-1/3 



(26) 



The observed width / D b s will be higher than this. First, because the observed 
width is the width of the projected synchrotron emitting rim, not the actual 
width. Second, since both advection and diffusion play a role, the synchrotron 
emitting shell has a width that is determined by the convolution of the ad- 
vection and diffusion length scales. So roughly l a b s ps ^/2l adv . The projection 
effects depend on the ratio of the physical width to the shock radius, and on 



possible deviations from spherical symmetry (e.g. Berezhko et al. 2003b I. 



Note that Eq. 26 is in principle less dependent on systematic uncertainties 
than either Eq. 20 and [23] because it does not depend on the measured 
velocity and it depends only weakly on the assumed compression ratio. But 
it does rely on the assumption that the emission is caused by electrons with 
energies close to the cut-off energy. This seems, in general, well justified for 
X-ray emission from young SNRs, because these are characterized by X-ray 
power-law indices that are steeper than the radio spectral index (a + 1). It 
is of interest to that Eq. [20] and [23] should only give similar magnetic field 
estim ates if the assumption that electr ons are close to the cut-off energy is 
valid (Vink 2005 Parizot et al. 2006). The fact that indeed the magnetic 



field estimate based on either Eq. 20 and 23 give similar answers ( |Ballet 



2006 ) gives additional justification for the use of Eq. 26 Interestingly, this 



also implies that the compression ratio cannot be deviating too much from 
r = 4, since Eq. 20 and 23 have a different dependence on r. However, if is 



there is a reason to expect that synchrotron radiation is observed at much 
lower frequencies than the cut-off frequency, Eq. [20] is to be preferred. See 
for example the use of Eq. [20] to det ermine the magnet i c field properties of 
a radio-arc in a cluster of galaxies by van Weeren et al. (2010). 
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Fig. 4 Overview of X-ray synchrotron emission profiles in several supernova rem- 
nants. The profiles are extracted from Chandra observations, with the exception 
of RX J1713. 7-39 46 which is based o n an XMM-Newton mosaic provided to us by 
Dr. Fabio Acero ( Acero et al. |2009b[ ). The solid red line indicates a best fit model 
based on a spherical model (i.e. trie emissivity is projected onto the sky) with an 
exponential fall off in emissivity away from the shock front. See text for further 
explanation. 
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2.6 Magnetic field measurements 



Fig. [4] shows the synchrotron emission profiles of several, relatively bright, 
supernova remnants. In all cases the emission is dominated by synchrotron 
emission, which means that for X-ray spectral line-rich supernova remnants, 
such as Cas A, Kepler, and Tycho, the line-poor spectral band of 4-6 keV 
was used. The red lines show emission models consisting of an emissivity that 
falls off exponentially toward the c enter (J oc exp(— (R — R s )/R w ), Rw is the 



width of the emitting region, c.f. Bamba et al. 2005), but, unlike Bamba 
et al. (2005), here spherical projection effects are taken into account. X-ray 



synchrotron emission from the pre-cursor is not fitted for here, because the 
lower magneti c field upstream of the shock should result in a much lower X- 
ray emissivity ( Berezhko et al. 2003b ). Nevertheless, this issue deserves more 



attention in future studies, as precursor emission can be used to constrain 
the diffusion constant upstream of the shock (see |Achterberg et al. 1994 for 
a study based on radio observations). 

As Fig. [4] shows, the model is far from perfect There can be several 
reasons why the observed profiles are not well fitted with the exponential 
model. First, in some cases there are large sca le deviations from sphe rical 
symmetry A case in point is SN 1006, for which Willingale et al. (1996) and 



Rothenflug et al. (2004) have shown that the X-ray synchrotron emission 
is coming from two caps, rather than from a spherical shell. This can be 
taken as evidence for a preference for acceleration in magnetic field orien- 
tations that are on average parallel to the shock normal ( |Rothenflug et al~ 
2004 ) . Indeed, Fig. [4] shows that the model is over-predicting the emission 



from the central region. Second, even though an exponential model may be 
a reasonable description in general, the acceleration of electrons to TeV en- 
ergies may fluctuate with time, and also the magnetic fields is likely to vary 
from one region to another. For steep synchrotron spectra (all spectra here 
have spectral indices of (r — 2 — 3.5) a s mall change in magnetic field lead 



to large changes in emissivity (I oc B 1 , Ginzburg and Syrovatskij 19671. 



Magnetic field variations can be the result of large scale Alfven waves, as 



explained in Sect 2.9 This caveat is related to the steady-state assumption 
that goes in deriving Eq. |20| For example, during a typical acceleration time, 
or synchrotron loss time, the shock velocity or the magnetic field may have 
substantially changed. Finally, magnetic-field amplification (see elsewhere in 
this volume and Sect. 2.7) operates most likely in the upstream region. In the 



downstream region the turbulent, amplified magnetic field may decay down- 
stream (Pohl et al. 20051. This would make in fact the magnetic estimates 



based on Eq. |20|an overestimate. But strong magnetic field decay would also 
affect the radio synchrotron morphology. This does not generally appear to 
be the case (e.g. |Cassam-Chena'i et al.| 2007), with the possi ble exception of 
one filament in the northeast of Tycho's supernova remnant ( Reynolds et al.| 
2011). Nevertheless, some form of decay may occur in general, and would 



affect the geometry of the X-ray synchrotron rims. 

4 Note that some regions were excluded from the fit, for example for Tycho (West) 
only the emission between 237" and 255" was fitted, as the emission closer to the 
center may also have a thermal bremsstrahlung origin. 
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Another assumption is that the spectral cut-off frequency is determined 
by radiative losses, i.e. we assume a loss-limited spectrum. These broader 
rims imply smaller magnetic fields (Eq 26 1 and, therefore, longer synchrotron 
loss times. During these synchrotron loss times the shock properties may have 
changed, i.e. the shock velocity was likely higher in the past, and the densities 
may have been different in the past. Another case where this assumption may 
not be valid is for the very young SNR Gl.9+0.3. 

We note here that a projected uniform thin shell emissivity fits in many 
cases the observed profiles as well, or even better than an exponential emis- 
sivity model. A thin shell has a peak that lies more inward of the shock 
radius. Indeed, the exponential model tends to underpredict the emission 
around the peak. A clear exception is Vela Jr (RX J0852. 0-4622), where an 
exponential model gives a much better fit than a uniform, projected shell. 
The best fit shell width of the uniform model is typically a factor two larger 
than the best fit characteristic width of the exponential model. Magnetic 
field estimates based on a uniform shell model are, therefore, 40% smaller. 

Table [2] lists magnetic field estimates of synchrotron rims in several super- 
nova remnants, based on Eq. 20 and the profiles shown in Fig. [4| In order to 
assess whether the spectra are loss-limited or age-limited we also list a rough 
estimate of the synchrotron loss times, assuming a typical photon energy of 
1 keV. It shows that only for Gl.9+0.3 (|Reynolds[ |2008| |Borkowski et al 



2010 Carlton et al. 2011) it seems likely that the spectrum is age- rather 



than loss limited. For this reason, one should treat the magnetic field esti- 
mates as upper limits. For the RCW 86, RX J1713.7-3946, RX J0852.0-3622 
and possibly SN1006 one should still be concerned about the steady-state 
assumption, as the synchrotron loss times are such that the shock velocity, 
magnetic fields may have changed during a time ~ t syn . 

The X-ray synchrotron spectra are caused by 10-100 TeV electrons that 
are also responsible for inverse Compton scattering background photons into 
TeV 7-rays. There is currently a controversy about whether the 7-rays emis- 
sion is dominated by these inverse Compton scattered background photons, 
or whether the 7-ray emission is due to the decay of neutral pions created 
by the collis ion of CR nuclei with background material (see 
for a review). The latter would directly estab 



maun 



2009 



Hinton and Hof- 



ish the presence 

of accelerated nuclei (hadronic CRs). In principle, once the average mag- 
netic field is known one can predict, for a given X-ray synchrotron flux, the 
expected inverse Compton TeV flux. 

An object that played a central role in the debate about the dominant 7- 
ray emission is RX J1713. 7-3946. It was by many considered to be an example 
of a supernova remnant emitting 7-rays dominated by pion decay ( Aharonian 



et al. 2004 Berezhko and Volk 



2008 



2010). But the lack of thermal X-ray 



emission could be taken as evidence for a low density inside the supernova 
remnant, which would argue against pion-decay dominated 7-ray emission 



and in favor of inverse Compton dominated emission (Katz and Waxman 



2008 Ellison et al. 20101. Another argument in favor of inverse Compton 



dominated 7-ray e mission is the close resemblance between the X-ray and 
7-ray morphology (Acero et al. 2009b I, which makes more sense if they are 



both caused by the same population of particles: electrons. Also the most 
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Fig. 5 

Dr. F. Acero (lAcero et al 



The XMM-Newton X-ray mosaic of RX J1713.7-3946 kindly provided by 
2009b I. The sector indicate (thin white line) indicates 



the region that was taken to make the emission profile in Fig. [4] 



recent Fermi-LAT results, combined with TeV results, favor inverse Compton 
dominated emission ( Abdo and Fermi LAT Collaboration]|2011| . This implies 
that the average magnetic field must be relatively small, 10 /iG. This seems 
at odds with the observed fluctuations of X-ray synchrotron filaments, which, 
if interpreted a s corresponding to syn chrotron loss time s, implies magnetic 

but see Sect. 



field of 0.1 mG (Uchiyama et al. 2007 



2.9 1. Often also the size 



of the X-ray synchrotron filaments are used to advocate a high magnetic field 
inside RX J1713. 7-3946. In that case one usually concentrates on bright X-ray 
synchrotron structures. But as the RX J1713. 7-3946 panel in Fig. [4] shows, 
these filaments are embedded inside a broad plateau of X-ray synchrotron 
emission. The overall X-ray synchrotron emission comes from a very broad 
region. For the south-eastern region of RX J1713. 47-3946 a magnetic field of 
B ~ 35 /iG is inferred (Table [2]) , but several other sectors of this supernova 
remnant have even broader emission zones (Fig. [5~|). An average magnetic 
field of B ~ 10 /iG is, therefore, quite plausible [ 



5 There has be en little attention to t his issue , but seejh ttp : //online . itp .ucsb . 
du/online/ astroplasmas_c09/vink/pdf /Vink_AstroPlasmasConf _KITP .pdf 
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Fig. 6 The downstream magnetic field pressure versus the energy flux through the 
shock poKs 3 - The data points are based on the values listed in Table [3] Note that 
the errors are quite large in both vertical and horizontal directions, and that the 
errors are largely systematic errors. The solid line indicat es the relation B 2 /(8ty) = 



10" 



poVf. The inset shows the values for SN1993J (Fransson and Bjornsson 



1998) and the solid line depicts the same relation as in t 



ae mam panel. 



2.7 The case for magnetic field amplification 



The average magnetic field in the interstellar medium is of the order of 5 pG. 
After a supernova remnant shock compresses the perpendicular magnetic field 
component with a factor four (or more) the resulting post-shock magnetic 
field should be B2 ~ 15 pG. As Table [2] shows, all supernova remnants seem 
to have magnetic fields higher than this; Cas A and Tycho have even con- 
siderably higher magnetic fields. This has been taken as evidence that some 
form of magnetic field amplification is operating in the vicinity of shocks of 



Bamba et al. 2005 Ballet 2006 I 



young supernova remnants (Vink and Laming 2003 Berezhko et al. 2003b 



Their are several means by which this magnetic field amplification occurs 
(see elsewhere in this volume), but a mechanism that received a lot of at- 
tention is the so-called Bell's instability (Bell 2004). For this mechanism the 
magnetic field energy density scales as B z oc u cr V s oc poVf, with u cr the CR 
energy density in the CR precursor. 



Volk et al. (2005) used their magnetic field determination based on the 
X-ray synchrotron rims to show that B 2 oc poV 2 , i.e. a roughly fixed fraction 
of the energy density in the shocked plasma comes from the magnetic field. 
However, Vink (2006), using a more up to date determination of the shock 



velocity of Cas A and slightly different magnetic field estimates, showed that 
B 2 oc poV? seems to fit somewhat better. 

Fig. [6] shows the downstream magnetic energy density versus poV^ for the 
magnetic fields listed in Tableland densities and velocities in Table [3] There 
is indeed a strong correlation. However, plotting B 2 versus poV 2 would give 
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a roughly similar correlation. The reason is that although there is a quite 
a dramatic contrast among young supernova remnants in p, there is little 
dynamic range in V s . Moreover, there is quite some systematic uncertainty 
in V s , as they have not been measured for all supernova remnants, and even 
if measured through proper motions, there is quite some uncertainty in the 
distance estimates toward the supernova remnants. What makes the case for 
a B 2 oc pV^ relation more plausible is that it better connects to the data 



point corresponding to the well studied radio supernova SN 1993 J (Fransson 



and Bjornsson |1998 Tatischcff 2009). This radio supernova was observed to 
have a very high shock velocity of 20,000 km s^ 1 , and its shock was moving 
through the dense wind of the progenitor. These values are dramatically dif- 
ferent from those of young supernova remnants. However, an assumption has 
to be that despite these different regimes for the parameters the underlying 
physics for the magnetic field amplification in SN1993J is the same as for the 
young supernova remnants. 

What is striking in Table [2] is that the lower magnetic fields are associated 
with larger supernova remnants (SN1006, RX J1713.7-3946, RX J0852.0- 
4622). This can be easily understood by noting that there is probably a 
linear relation between B 2 and po and that for X-ray synchrotron radiation 
one needs shock velocities > 2000 km s^ 1 . Only supernova remnants in 
low density environments can reach radii of more than 7 pc and still have 
such high velocities. But low densities does, according to our knowledge of 
magnetic amplification, correspond to low magnetic fields. 



2.8 X-ray synchrotron flux changes 



As a supernova remnant expands its shock velocity decelerates and the 
plasma and relativistic particles cool adiabatically. The synchrotron emis- 
sion will, therefore, change as a function of time. This was first noticed for 
the radio-synchrotron emission from Cas A whose flux density changes by 
about 1% per year (Hogbom and Shakeshaft 1961 Baars et al. 1977). This 



secular decrease was interpreted by Shklovsky (19681 as due to adiabatic 
losses. 

For X-ray synchrotron radiation adiabatic losses may also be taken into 
account, but a more prominent factor is the change in shock velocity, as this 



changes the cut-off frequency in the loss limited case (Eq. 19), which has a 
strong effect on the synchrotron flux above the cut-off frequency. If we take 
for example the following asymptotic form of the synchrotron spectrum (c.f. 
Zirakash vili and Aharonian] |2007| |Katsuda et al.| |2010b| |Patnaude et al.[ 
20111): 



Fihv) = <j)B r ° x (hv) 



exp 



t'cut-off 



with Po the asymptotic spectral index below the cut-off and < 
constant, we see that the spectral index is given by: 



P 




(27) 

> a normalization 
(28) 
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Fig. 7 The flux decrease of the SNR Cas A in the 4.2-6 keV, synchrotron - 



dominated, X-ray band, based on Chandra observations (Patnaude et al. 
(Courtesy Dr. D. Patnaude, reproduced by permission of the A AS. 



2011) 



with fa'cut-off given by Eq. 19 



The fractional change in cut-off frequency as a function of the change in 
shock velocity is determined by: 



1 ^cut-off _ - o ( m ~ 

t'cut-off dt V s dt t 



(29) 



with t the age of the supernova remnant and m the so-called expansion pa- 
rameter characterizing the evolution of the supernova remnant (R oc t m , V s — 
mR/t). 

We can use this to estimate the expected X-ray synchrotron flux change 
due to changes in hi/ cu t- s alone|^J 



1 dF{v) 
F{v) dt 



1 



dv. 



cut-off 



cut— off 



dt 



"(m-1) 



t'cut-off 



(30) 



A recent analysis of the Chandra observations of Cas A showed that 
this supernova re mnant has a decl i ne in X-ray synchrotron flux at a rate 
of 1.5%-2 % yr- 1 (|Patnaude et al.||2011 . see Fig. ff\. In fact, if one takes 



m = 0.65 ( |Vink et al.[|1998[|Delaney and Rudnick 2003 Patnaude and Fesen 



Note that 



Katsuda et al 



(2010b I uses a similar formalism, but with slightly 
different dependencies as they also consider the change in the amplified magnetic 
field 
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2009), t — 330 yr, and ■J v/vcut-oS ~ 0.6 — 7 (corresponding to /if cu t-off 



0.1 — 2 keV), then Eq. 30 predicts a fractional decline rate of 0.06-0.7% yr 



This is smaller than observed, but it should be noted that the strongest 
decline rate seems to be associated with the X-ray synchrotron emission 
from the reverse shock (Sect |2~1~2 | for which m is not well measured. We 
also note that the cut-off frequency does not seem to be well defined. Indeed, 
one of the peculi ar features of Cas A it s near power-law spectral shape from 

whereas Eq 27 predicts a gradual 



2006) 



- 5 - 100 keV ( jRenaud et al. 
steepening. 

Eq. [28] can be used to couple the flux decline rate with a steepening of 



the synchrotron spectrum (Patnaude et al. 2011): 



1 dF{v) 
F{v) dt 



= -2 



dt' 



(31) 



The Chandra observations of Cas A indicate dFjdt = 0.022 yr -1 , which 
should correspond with a flux decline of ~ 4 % yr -1 , which is clearly too 
large. These discrepancies could either point to the invalidity of the simple 
model for flux decline (see also the discussion in Patnaude et al. 2011), or 



the invalidity of Eq [27] to describe the overall X-ray synchrotron spectrum of 
Cas A. 

For the youngest known supernova remnant Gl.9+0.3 the X-ray syn- 
chrotron flux is measured to be increasing with a rate of 1.7 ± 1.0 % yr -1 
( |Carlton et al. 2011). A similar flux increase is observed in the radio. This 
flux increase is most likely due to an expansion in a uniform medium. Due 
to an increase of the swept up mass the emission will increase, although not 



linearly with the swept- up mass due to adiabatic losses (Carlton et al. , 2011 ). 



The same is true for the X-ray synchrotron emission, but here, like for Cas 
A, the change in cut-off frequency will be a large effect. If the synchrotron 
spectrum is age limited the change in cut-off frequency can be found by 
combining Eq. 
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and Eq. 18 which gives (Katsuda et al. 2010b): 



^cut 



_ off cx B 3 V S H 2 , 



which, if B is assumed constant, gives: 
1 cut-off 4(to-1) + 2 



^cut-off 



dt 



t 



(32) 



(33) 



This implies that the cut-off frequency will decrease for m < 0.5, and con- 
sequently the X-ray synchrotron flux would be more likely to decrease. Note 



of B 2 cx V? and find 



that Katsuda et al. (2010b) and Carlton et al. (2011) include a dependence 



1 



dp, 



cut— off 



(7(m-l) + 2) 



^cut-off dt t 

with a decline in cut-off frequency for m < 0.71. For B 2 oc K, 3 we find 
1 dfcut-off (17(m- l)/2 + 2) 



^cut-off 



dt 



t 



(34) 



(35) 
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which gives a decline in cut-off frequency for to < 0.76. These expressions can 
be used together with the second part of Eq. 30 to estimate flux changes in age 
limited cases. Note that the value for to is not known for Gl.9+0.3, because 
it is usually estimated from the relation V s — mR/t, but since Gl.9+0.3 is 
not connected to an historical event, t is not known. 



The equations 33 - 35 show a peculiarity of the X-ray synchrotron emission 



for age limited spectra. If the m is too small, the cut-off frequency only 
decreases. For the constant magnetic field case this is not so dramatic because 
young supernova remnants in Table [3] have in most cases m > 0.5. But in 
the case of magnetic- field amplification it may in some cases be difficult to 
understand why there is X-ray synchrotron radiation at all. For example, for 
Cas A, to = 0.65. This value corresponds roughly with a Sedov model for 
a supernova remnant evolving inside the wind of the progenitor. Since Cas 
A evolves in a dense w ind to — 0.65 is probably valid for the last 200 yr 
( van Veelen et akj 2009 Fig. 7). If that is indeed the case then most of the 
electrons responsible for X-ray synchrotron radiation may have been electrons 
that were injected in the early supernova remnant phase. Indeed, in the early 
phase one expects to > 0.8. Fo r example, SN1993J was observed to have 
to 0.84 (Marcaide et al. 



2009). 



For Type la supernova remnants a similar problem exists if one takes 
them to evolve in their early phase in a self-similar way. One would then 
expect to = 0.57 (Chevalier 1982a). But it seems more likely that the eje cta 
structure has an exponential profile (e.g. Dwarkadas and Chevalier 1998), in 
which case in the early phase one expects to > 0.7. 

One should note, however, that magnetic-field amplification is not an 
instantaneous process. It probably evolves simultaneously with the acceler- 
ation of particles to higher and higher energies. This, nevertheless, shows 
that observations of X-ray synchrotron flux decline serve to cast a new light 
on the theory of the evolution of supernova remnants, the theory of X-ray 
synchrotron radiation, and even on the issue of magnetic field amplification. 



2.9 X-ray surface brightness variations in supernova shells 

Clumpy structures in X-ray emission are frequently observed in imaging ob- 
servations of supernova remnants. These structures have X-ray spectra that 
can be both thermal or non-thermal in nature. Structures with thermal X- 
ray spectra may be caused either by instabilities at the contact discontinuity, 
by clumps of ejecta protruding beyond the blast wave or by reverse shock 
heating of clumpy ejecta material. The spatial and temporal characteristics 
of structures with non-thermal X-ray spectra are probably caused by par- 
ticle acceleration and the associated processes of fluctuating magnetic field 
amplification. 

Here, we limit ourselves to a discussion of X-ray brightness variations 
in supernova shells related to energetic particle acceleration processes and 
therefore to the structures with non-thermal X-ray spectra. These X-ray 
brightness variations may originate both from the bremsstrahlung emission of 
electrons with keV-MeV energies and from synchrotron radiation of electrons 
with TeV energies in stochastic magnetic fields. 
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2.10 Synchrotron surface brightness fluctuations and stripes in Tycho 



X-ray synchrotron emission structures have been observed with the superb 
spatial resolution of the Chandra telescope in many young supernova rem- 



nants (see e.g., Vink and Laming] |2003| |Bamba et al.||2005||Uchiyama et al. 



2007||Patnaude and Fesen||2009[|Helder and Vink[ 2008; E riksen et al.[|2011[ ). 



The morphology of the extended, nonthermal, thin filaments observed at the 
supernova remnant edges, and their X-ray brightness profiles, strongly sup- 
port the interpretation that > 10 TeV electrons are accelerated at the 
forward shock of the expanding supernova shell and produce synchrotron ra- 



Reynolds 


2008 


Vink 


2008a 



The small scale X-ray structures brightening and decaying on a one-year 



pernova remnant RX J1713. 7-3946 and in Cas A by (Patnaude and Fesen 


2007 


Uchiyama and Aharonian 


2008 


Patnaude and Fesen 


2009). The rapid 



variability of the structures detected with Chandra were interpreted as the 
X-ray synchrotron structures that are produced by ultrarelativistic electrons 
accelerated in a strongly magnetized environment with the magnetic field 
amplification by a factor of more than 100. Moreover, Uchiyama et al. (2007) 



suggested that the X-ray variability is a direct signature of the ongoing 
shock-acceleration of electrons in real time. However, the presence of the 
very high mag netic fields i n the shell of supernova remnant RX J1713. 7-3946 
was shown by Butt et al. ( 2008 1 to be too constraining from the view of the 



multi-wavelength observations, unless the high field regions occupy a small 
fraction of the volume. 

Efficient diffusive shock acceleration of high-energy particles requires a 
substantial amplification of magnetic field fluctuations in the vicinity of the 
shock; see e.g. |Bell| p978[); |Blandford and Eichler| (fl987|); |Malkov and Drury 



(2001). Magnetic-field amplification mechanisms due to CR instabilities in 
nonlinear diffusive shock acceleration were proposed recently by several stud- 
" "2005]); |Amato and Blasi] ([20061; IvTadimirov et "aL[ ([20061; 



ICS 



Bell (2004 



Marcowith et a 



■| H2006D; IVladimirov et al.] Q2009D; |Zirakashvili and Ptuskin 



Bykov et al.| ( |2011b| ); |Schure andBell| ( |201ip . These models predict 



p008[ ) _ 

cosmic-ray streaming ahead of the shock will amplify the local magnetic field 
substantially. 

An alternative explanation for the surface brightness fluctuations in RX 
J1713. 7-3946 and Cassiopeia A was offered by |Bykov et all (|2008b|). As the 



amplified magnetic fields in supernova shells arc highly turbulent, they might 
produce local enhancements, resulting in isolated structures in a synchrotron 
image that show large variations in the brightness. Note that these brightness 
variations can occur even if the particle distribution is smooth and steady. 
Bykov et al. ( 2008b| 2009) investigated this effect, resulting in the follow- 



ing three predictions. Firstly, the cut-off energy of the emitted synchrotron 
spectrum can be significantly increased for a distribution of electrons residing 
in a turbulent magnetic field, as compared to electrons in a uniform field with 
the same strength. Secondly, the variation in surface brightness is on shorter 
timescales for images obtained at higher energies. More specifically, the time 
scales expected are similar to those found in X-ray images of RX J1713.7- 
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3946 and Cassiopeia A. Future hard X-ray imaging instruments as NuSTAR 
and ASTRO-H will probably be able to test this prediction. Thirdly, the 
model predicts that the synchrotron emission from these bright structures 
will be highly polarized, as illustrated in Figure [HJ 

Recently, very unusual structures consisting of ordered sets of bright, 
non-thermal stripes were discovered a deep Chandra exposure of Tycho's 
supernova remnant (Eriksen et al. |2011 ). The stripes are clearly seen in the 
4.0-6.0 keV image of the remnant (Fig. [9]). At these energies, the emission is 
dominated by X-ray synchrotron radiation. Understanding these structures 
presents a unique opportunity for current models of X-ray synchrotron images 
of young supernova remnants. As the stripes are likely caused by peaks in 
magnetic turbulence of a perpendicular shock, this provides a way to measure 



the orientation of the magnetic field of the ambient medium ( |Bykov et al. 



2011a 



Also, the coherent nature of the X-ray stripes likely suggests that the 
underlying magnetic turbulence is strongly anisotropic relative to an or- 



dered mean background magnetic field. |Bykov et al 



cho's stripes assuming they are the result of the non 



(2011a) modeled Ty- 



inear evolution of the 



anisotropic CR-driven magnetic instability by ( Bell 2004 Vladimirov et al 



2009 Bykov et al. 2011b I, and concluded that the structures should only 



occur in the shell section where the local field lies mainly along the shock 
surface and where the turbulence cascading is suppressed. The simulated 
stripes are shown in Fig. 10 They found that the orientation of the ambient 
magnetic field at the location of the stripes has to be quasi-perpendicular 
to the shock normal. The unstable growing magnetic modes must maintain 
coherence over a scale, close to the size of the CR precursor. Additionally, 
the modes have to be linearly polarized. Turbulence cascading along the 
mean field should be suppressed to prevent the broadening in wavenumber 
(k) of the generated turbulence. Vladimirov et al. (2009) demonstrated in 
their nonlinear diffusive shock acceleration model that the spectral cascade 
suppression results in a peaked structure of the magnetic fluctuation spectra. 
The predicted polarized fraction of ~ 50% would make these stripes easily 
detectable in future X-ray polarization observations 



2.11 The asymptotic behavior of particle spectra 

According to non-linear diffusive shock acceleration theory the spectrum of 
accelerated particles should not be a power law, but instead gradually harden. 
The reason is that the presence of accelerated particle ahead of the shock 
creates a shock precursor, which pre-compresses the gas, before it is further 
compressed by the gas shock. Also the overall shock compression can be 
substantially higher than the canonical factor four, for a high Mach number 
shock in a monatomic gas. At the same time, a high compression in the pre- 
cursor, means a lower Mach number at the actual gas shock, as the medium 
has already been adiabatically heated. All this means that the highest energy 
particles may scatter back and forth across the total shock region, sampling 
the full compression ratio, whereas the low energy cosmic rays only sample 
the gas-shock region, which may have a compression ratio of about r — 2.5. 
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Fig. 8 Simulations of polarized synchr otron emission in a random magnetic field at 



5 keV adapted from Bykov et al. ( 2009 1. The left panel shows intensity. The central 



panel shows the product ot intensity and polarization degree. The right panel shows 
the degree of polarization (colorbar represents its scale). The \J {B 2 } — 3 x 10~ 5 
G and the spectral energy density of mag netic fluctuations W(k) oc k~ 2 . For the 
models with W{k) oc AT 1 see Figs. 2-4 in Bykov et al.| |2009|). 



Since in the test particle approach the power-law index of the spectrum of 
accelerated particle is q = (r + 2)/(r — 1) one sees that for r = 2.5 q = 3, 
whereas for r = 7 it is q = 1.5. For non-linear shock acceleration one expects 
a spectrum that gradually changes in spectral slope as a function of energy. 
For low energies the slope of the particle spectrum will be q w 3, but it is 
expected to asym ptotically harde n to q = 1.5 at high energies if compression 
ratios are r « 7. Malkov ( 1997 ) showed that for even higher compression 
ratios one expects q = 1.5. On the other hand, if Alfven waves are on aver- 
age directed away from the sho ck front a softer spectrum may be expected 
( Zirakashvili and Ptuskin 2008 1 . 



Observationally, there is clear evidence that synchrotron spectra from 
young SNRs are not power laws over a long range of frequencies. The first 
evidence for this was presented by Reynolds and Ellison (1992), who showed 
that the radio spectra of Tycho and Kepler gradually harden from ~ 10 8 
to - 10 10 Hz. 



Hz 



Cas A has a rather steep radio spectrum, a — 0.78, corresponding to 
q = 2.6 and, on face value, corresponding to a compression ratio of r = 2.9, 
with little evidence for intrinsic spectral hardening within the radio spectrum. 
But comparing the radio data with synchrotron emitting regions detected 
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in the infrared reveals that the synchrotron spectrum must have hardened 



(Jones et al. 2003 1 . 



Also some X-ray synchrotron spectra show that spectral hardening must 
occur. In X-rays the synchrotron emission is observed near or beyond the 



spectral cut-off frequency, but for RCW 86 (Vink et al. 20061 and SN 1006 



( Allen et al. 2008 ) the extrapolation of the radio synchrotron spectrum can- 



not at the same time explain the X-ray flux and the X-ray spectral index. For 
both cases a curved synchrotron spectrum describes the combined radio and 
X-ray observations better. Note that in both cases an heuristic curved elec- 
tron spectrum was chosen. In Fig.[TT]thc data points and curved synchrotron 
models are shown for the northeastern region of RCW 86. Both models with 
an asymptotic power-law index for the electron distribution of q = 2 (the 
classic Fermi acceleration index) and q — 1.5 can in principle explain the 
observations. Note that incorporating curvature in the models will lead to 
lower values for the cut-off frequency when compared to a classic, but very 



much simplified, X-ray synchrotron model as srcut (Reynolds and Keohane 
1999). Note that one should take into account that the radio emission from a 
SNR is likely to come from a much larger volume than the X-ray synchrotron 
emission. So it is better to apply this type of analysis to small X-ray emitting 
regions, rather than to emission from the entire SNR. 

Although the evidence for curvature is quite compelling, we see that the 
evidence for an asymptotic particle power-law index of q — 1.5, expected for 
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Supernova remnant synchrotron emission images simulated by Bykov et al. 

Simulations done within 



Fig- 10 

(2011a I to model the Tycho supernova remnant stripes 



a trame of the nonlinear model of diffusive shock acceleration taking into account 
magnetic field amplification from a CR current driven instability. The left panel is 
the synchrotron X-ray intensity at 5 keV regime. The degree of polarization of the 
X-ray emission is in the right panel with the degree of linear polarization shown in 
the color bar. The high degree of polarization of the X-ray synchrotron stripes is 
evident. Reproduced by permission of the AAS. 



strongly non-linear acceleration, is much less clear. Also TeV 7-ray obser- 
vations do not provide evidence for spectra that are that hard, with most 
observed spectral indices r > 2. But the measurements errors are still large 
enough, and the interpretation of the results (pion-decay or inverse Comp- 
ton scattering) still uncertain enough, that one cannot yet firmly rule out 
q = 1.5. 
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Fig. 11 Broad band energy density spectrum of the northeastern part of RCW 
86. The radio spectrum, here represented with a flux density point at 1 GHz and a 
spectral index a = 0.6, cannot be connected with the X-ray flux point if one assumes 
an electron spectrum consisting of a power law with an exponential cut-off. Instead 
a gradual flattening is needed, although there is enough freedom to either have a 
modest flattening, with an asymptotic index of q — 2 or to have q = 1 -5 predicted 



by non-linear shock acceleration theory. (Figure from Vink et al. 
by permission of the AAS 



2006 reproduced 



2.12 Particle acceleration at the reverse shock 



Throughout the Sedov phase of the supernova remnant evolution, the shock 
velocity of the reverse shock can be comparable to, or even higher than, the 



forward shock velocity (Truelove and McKee 1999). This in itself makes the 



reverse shock a probable location for CR acceleration. However, in contrast 
to the magnetic field of the ambient medium, the magnetic field in the freely 
expanding ejecta (which is upstream for the reverse shock) is negligible. As 
magnetic fields are essential for particle acceleration processes, it is not obvi- 



ous that reverse shocks are able to accelerate particles (Ellison et al.l 2005). 



Observational evidence for cosmic-ray acceleration at the reverse shock would 
provide useful constrains on the required magnetic field upstream to make 
particle acceleration work. 

A number of observations suggest synchrotron emission from the reverse 
shock. DeLaney et al. ( |2002[ ) found radio synchrotron emission coinciding 
with the reverse shock, implying the presence of GeV electrons. Also, the 
southwest corner of the RCW 86 remnant shows signs of X-ray synchrotron 
emission. This corner is interacting with a cavity wall and as a reaction, the 
reverse shock probably developed very strongly in this part of the remnant. 
Coincident ally, it is exactly in this part of the remnant where there are dis- 
tinct filaments of non-thermal X-ray emission inside of the thermal emission 
at the outside ( Rho et al. 2002 ) , suggest ing the presence of accelerated par- 
ticles at the reverse shock. Additionally, |Zirakashvili and Aharonian ( 2010 1 
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suggested that the arc-like structure in the X-ray synchrotron emission of 
RX J1713. 7-3946 is formed by the reverse shock as well (See Figure [5]). 

The non-thermal X-ray emission of Cas A is rather complex. The ori- 
gin of the narrow filaments that mark the outside of the remnant is at the 



forward shock, bu t there is non-thermal em ission at the inside as well ( Got 
thelf et al. 2001). Helder and Vink| (20081 concluded from a deprojection 



of the megasecond Chandra image in the 4 to 6 keV band and the different 
kinematical properties ( DeLaney et al. , 2004[ ) that these inner non-thermal 
X-ray filaments are originating from a shell inside of the forward shock, rather 
than located at the forward shock , projected onto the center of the remnant. 
Uchiyama and Aharonian (2008) reached a similar conclusion for Cas A, 



based on the fact that the inner filaments showe d surface brightness fiuctua - 
tions and the outside structure did not (but, see Patnaude and Fesen 2009). 
A non-thermal bremsstrahlung is unlikely for these inner filaments, as non- 
thermal tails to the electron distribution would have disappeared through 
Coulomb interactions on timescales shorter than the lifetime of the remnant 
dVinkl |2008bl. Therefore, iHelder and Vinkl (|2008b concluded that the non- 



i 1 1 i - 1 i i i 

thermal X-ray emission at the center of Cas A is synchrotron radiation of 

TeV electrons at the reverse shock. Whether these electrons were accelerated 
or reaccelerated (ISchure et all 120101) at the reverse shock remains unclear. 



2.13 Gamma-ray observations of supernova remnants 



Ever since the detection of 7-rays from the Crab nebula by the Whipple Tele- 
scope ( Weekes et al. 1989 ) showed the potential of Cherenkov telescopes, the 
field of TeV 7-ray astronomy is thought to hold the key to solve the mystery 
of the origin of Galactic CRs. But despite the tremendous advances in TeV 
7-ray astronomy, there is, as we will discuss below, no definitive proof that 
the CRs are accelerated by supernova remnants in sufficient quantities and 
to energies beyond the knee in order to explain the CR spectrum on earth. 
Nevertheless, our understanding of CR acceleration in supernova remnants 
has increased substantially as a result of 7-ray astronomy. The results of TeV 
7-ray telescopes (~ 0.1 — 100 TeV) has recently been complemented by ob- 
servations w ith the GeV 7-ray satellites (~ .1 - 100 GeV) AGILE ( jTavani 



et al.. 2008) and Fermi (Atwood et al. 2009). 



2.13.1 Hadronic versus leptonic ^-ray radiation 



7-ray astronomy is of great importance for determining the CR content of 
astrophysical sources, because only in 7-rays one can observe photons emitted 
as a result of hadronic CRs (i.e. CR nuclei/ions), which make up 99% of 
the CRs observed on earth. The emission is caused by CR nuclei colliding 
with the background medium, thereby producing, among others, pions (pi- 
mesons). The neutral pions decay immediately into two photons, each having 
in the rest-frame of the collision an energy of half that of the pion rest mass, 
M^oc 2 = 135.0 MeV. In the frame of the observer the photon has an energy 
that is on average about 12% of the energy of the primary CR particle. The 
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7-ray spectrum and emissivity is directly related to the spectrum of CRs in 
the source and the density of the background medium, tih- If the primary 
particle spectrum is given by 



n CI (E p )dE p = kE p v *dE p 



(36) 



with k a normalization constant and q p the spectral index, the 7-ray emis- 
sivity is given by 



dn(hv) 
dhvdt 



= n a 



ck n E, 



-qp 



pp,incl 



(E p ) 



dfjhv, Ep) 
dhv 



dE„ 



(37) 



Here c is the velocity of the CR, <7 PP} i nc i the inelastic collision cross sec- 
tion, and the function df(hv,E p )dE p /dhi> gives the energy distribution of 
the 7-rays as a function of the primary proton energy. This function includes 
the effects of multiplicity (the number of pions produced per collision) , the 
distribution of energies of the pions, and emission angle dependence. In ad- 
dition, it should include the effects of composition of both the cosmic rays 
and the background plasma. An approximation that is sometimes made is 
that for a large range of energies <7 PPi i ne i is assumed to be constant , and 
< hv/E p >» 0.17 ( including the contribution of 77-meson production, Aha 
ronian and Atoyan] 2000). S ince the pion prod uction cross section is only 

for approximate formulae), 



mildly energy dependent (see Kelner et al. 2006 



the spectral photon index is approximately that of the spectral index of the 
particles, r s» q p . 

However, the detection of 7-rays cannot be uniquely attributed to CR nu- 
clei, because relativistic electrons (leptonic CRs) also produce 7-rays through 
two radiative processes: inverse Compton scattering and bremsstrahlung. 

The emissivity of inverse Compton scattering depends on the density of 
relativistic electrons and the number density and spectrum of the background 
photons that are being upscattered. The relation between initial photon en- 
ergy and up-scattered energy is hi^c ~ lt^ v i- The background radiation is 
often assumed to be from the cosmic microwave background, which usually 
dominates the photon field in the Galaxy. But in some cases other sources of 
photons dominate. A case in point is Cas A, which is itself a strong source 
of far infrared radiation, which is the dominant source of seed photons for 



inverse Compton 7-rays in the object (e.g. Abdo et al. 2010b ) 



The cross section of inverse Compton scattering is given by the Thomson 
cross section ctt, except at very high energies where the cross is reduced as 
a result of the Klein-Nishina effect. Thus the total electron energy loss due 
to inverse Compton scattering is approximately 



dE c 
dt 



-in 



-0Tit rad 7 o 



(38) 



with u ra d the radiation density of background photons. Since vie oc E% (so 
dE oc v~ 1 / 2 di') 1 the relation between the photon index and the electron spec- 
tral index is r = (q c + l)/2, with q e the spectral index of the electrons. This 



7 Sometimes the effect of multiplicity and cross section is comb ined in an inclusive 
cross section a(i? p ,i n ci) (e.g. Dermer 19861 iKamae et al. 20061. 
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Fig. 12 7-ray excess count map of RX J1713. 7-3946 as observed by HESS. Over- 
layed are 1-3 keV X-ray (ASCA) intensity contours. Credit: Aharonian et al., A&A, 
235, 243, 2007, reproduced with permission © ESQ. 



relation between the spectral indices of the photon and electron spectrum 



is the same as for synchrotron radiation ( jGinzburg and Syrovatskii 19651, 



and is quite distinct from the relation between primary spectrum and photon 
spectrum expected from pion decay. Note that over a large range the elec- 
tron and proton spectral indices are expected to be similar with q e = q p m 2. 
The inverse Compton spectrum is expected to be harder F « 1.5 — 1.7, but 
in the TeV range of the spectrum the cut-off of the electron spectrum can 
still result in a spectrum softer than r = 1.7. For that reason, the broad 
coverage provided by including GeV observations with Fermi or AGILE is 
important. Moreover, the GeV emission of supernova remnants comes from 
the same electrons that also produce radio synchrotron radiation. So the 
spectral slopes in the radio and 7-rays can be directly compared. 

A third radiation process that is important, and which like pion-decay 
scales with background density nn, is bremsstrahlung, caused by deflections 
of the electrons as they encounter charged particles. For relativistic electrons 
the power-law spectral photon index is approximately equal to the primary 
electron spectrum, r — q c . The bremsstrahlung is approximately a power 
law over a broad range of energies, unlike pion decay, which drops off rapidly 
for energies below half the pion rest mass energy of neutral pions 65.5 McV. 

2.13.2 j-ray observation of young shell-type supernova remnants 

The first shell-type supernova remnant to be detected in ultra high energy 



7-rays was Cas A with the HEGRA single-dish Cherenkov telescope (Aha 
ronian et al.| [2001). With the advent of more sensitive telescopes like HESS 



and later MAGIC and VERITAS, more shell- type supernova remnants were 
detected (Table |1.1| . These include well known sources such as Cas A ( Aha- 



ronian et al.| [200l| | Albert et aH[2007b| , SN1006 ( |Acero et al.| |20Tol) RCW 
86 (Aharonian et al. 2009) and very recently Tycho ( |Acciari et al.| 20111. 
In addition, there are a number of previously poorly known supernova rem- 
nants detect ed in TeV 7-rays. Thes e include RX J1713. 7-3946 (Aharonian 
eTaLl [20041 and RX J0852.0-4622 1HESS Collaboration: F. A. Aharonianl 



2006). These two sources have a low surface brightness in the radio and X- 



rays, but have turned out to be among the brightest TeV sources in the 
Galaxy. All these TeV sources are also known as X-ray synchrotron emit- 
ters, with the X-ray emission from RX J1713.7-3946 and RX J0852.0-4622 



being even completely dominated by synchrotron emission (Sect. 2.3). The 



average power-law index of the 7-ray emission from young shell-type super- 
nova remnants is _T 2.3 ± 0.3. The hardest emission was detected for RX 
J1713. 7-3946, with emission detected up to ~ 100 TeV. Its average spectrum 
is best fit with r = 2.32 ± 0.01. But a better fit is provided by a power law 
with an exponential cut-off that scales as exp(— \J v/v c ). In that case the best 
fit parameters are r = 1.79 ± 0.06, a cut-off energy of hv c = 3.7 ± 1.0 TeV. 

Initially a strong case was made that the TeV 7-ray emissio n was caused 
by pion-decay radiation (e.g. Berezhko and Volk 2008, 2010). This would 



constitute direct proof that proto ns a re accelerated to energies of at least 
~ 10 14 eV. As discussed in Sect. 2.3 this implied that the magnetic field 
must be relatively high in all these sources (> 100 /iG), as otherwise the 
observed synchrotron flux would require a high CR electron density. This 
would enhance the 7-ray inverse Compton scattering, and making it the 
dominant source of 7-rays. 

The fact that all young supernova remnants that emit TeV 7-rays are 
also X-ray synchrotron emitters (Table is a reason to be cautious about 
pion-decay dominated models, since, clearly, electrons are present with en- 
ergies > TeV, which could result in an appreciable, if not dominant, inverse 
Compton contribution to the 7-ray flux. Moreover, the lack of thermal X-ray 
emission from the prominent 7-ray source RX J1713. 7-3946 suggests that 
the density in the supernova remnant is low. This would suppress pion-decay 
emission, as it scales with nn (Eq. |37j see Katz and Waxman 2008 ; Ellison 
2010 1. But it may be that the thermal X-ray emission is suppressed 



et al. 



by very low plasma temper atures, which coul d be a by-product of e fficie nt, 
non-linear CR acceleration (Drury et al. 2009). As explained in Sect. 2.3 the 



magnetic field in RX J1713. 7-3946 may not be as high as sometimes assumed. 
This makes it also more likely that inverse Compton scattering dominates 
the TeV emission from this source. 

The best evidence that the TeV 7-ray emission from RX J1713. 7-3946 
is dominated by inverse Compton scattering is provided by recent Fermi 
observations (Abdo and Fermi LAT Collaboration 2011), because they pro- 



vided the broad spectral coverage needed to distinguish the different spectral 
behavior of pion-decay and inverse Compton scattering. This showed that 
r = 1.5, which is consistent with the spectral index expected for inverse 
Compton scattering. The hadronic-model (pion-decay) for the 7-ray is not 
yet completely ruled out, as r — 1.5 is still consistent with a very hard pri- 
mary proton spectrum of q p = 1.5, which could be the result of extremely 
efficient, non-linear CR acceleration (Malkov 1997). But in that case the 
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spectral break around 4 TeV indicates that the maximum proton energy is 
around 40 TeV, very far removed from "the knee" . So ironically, a leptonic 
model may be more consistent with the idea that supernova remnants accel- 
erate protons to very high energies, as it gives the freedom to assume that 
the non-detected protons have much higher cut-off energies. 

Cas A was the first shell-type supernova remnant detected in TeV 7-rays 



(Aharonianet al. 2001| Albert et al.[ 2007b). Its broad band GeV/TeV 7-ray 



spectrum casts some doubt on the idea that more than 10% of the explo- 



sion energy in young supernova remnants is contained by CRs (Abdo et al. 



2010b). The recent Fermi observations do not provide sufficient information 



on the dominant 7-ray emission, as both hadronic (pion decay) models and 
leptonic models (for Cas A a combination of inverse Compton scattering and 
bremsstrahlung) provide reasonable fits to the 0.1 GeV-10 TeV spectrum. 
But for Cas A the conclusion is that either model shows that the CR energy 
budget in Cas A is at most 4 x 10 49 erg, which corresponds to at best 2% of 
the total explosion energy. 

Of all the young supernova remnants, only Tycho's supernova remnant 
seems to offer a good case that a substantial amount of the explosion energy 
has been used to accelerate CRs, and that the 7-ray emission is dominated by 
pion-decay. Tycho's supernova remnant was only recently detected in the TeV 
and GeV range by res pectively VERITAS (Acciari et al. 2011) and Fermi 
(Giordano et al. 2011 ). The spectrum in the GeV range has a spectral index 
of r = 2.2, whereas the radio spectral index a = 0.6 implies an inverse 
Compton scattering spectrum of -Tic — 1.6. Non-thermal bremsstrahlung 
could in principle fit the spectrum, but requires a density of rtn ~ 9 cm~ 3 . 
This would be much higher than expected based on the dynamics of the 
supernova remnant (Katsuda et al. 2010a), unless the compression ratio 



behind the shock front is very high. The total energy in CRs required by 
fitting a pion-decay model to the joint Fermi and VERITAS data shows that 
10 ± 5% of the explosion energy is contained in CRs. The uncertainty is 
caused by some uncertainties in the distance, density, and overall explosion 
energy. 

So what are we to make of the 7-ray emission from young shell-type 
supernova remnants? Taking all the 7-ray observations of young supernova 
remnants together, the results are somewhat mixed and confusing. Sources 
that initially were considered promising candidates for substantial hadronic 
CR acceleration, RX J1713. 7-3946 and Cas A, do not seem to require a 
substantial energy in CRs. But Tycho's supernova remnant seems to have 
a substantial energy in CRs, although perhaps not as high as one might 
hope for in a young supernova remnants. Currently, only a few TeV emitting 
supernova remnants have also been detected by Fermi, so the balance may 
still shift a little more toward a hadronic origin for 7-ray emission. 

However, an important ingredient of CR acceleration by supernova rem- 
nants has only recently received more attention. This is the issue of CR 
escape. Non-linear CR acceleration is thought to be accompanied by sub- 



stantial escape of high energy particles ups tream (e.g. |Malkov and Drury 
[200Tj |Blasi et al.| |2(J55l [keville et al.[ |200^1 |Vink et al.j |2010| |Drury[ |2011 



the latter for a more nuanced view of escape). It is actually this escape of CRs 
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that makes the high compression ratios possible. It is therefore possible that 
for sources like Cas A a substantial amount of CR energy has already escaped 
the supernova remnant. Also the discovery of TeV sources in the vicinity of 
mature supernova remnants testifies of the importance of CR escape (see the 
next section). 

Note that core-collapse supernova remnants such as Cas A and Type 
la supernova remnants, such as Tycho, show a different 7-ray response to 
escaping CRs. The reason is that core-collapse supernova remnants often 
evolve in the wind of their p rogenitors, which in the ca se of Cas A is likely 
a dense red supergiant wind ( Chevalier and Oishi{ 2003[ ). In such a wind the 



density falls off with radius as p(R) — M w /(4ttR v w /3), or, nn(R) oc R~ 
So escaping CRs, while diffusing outward, encounter less dense material. 
If we assume that the total number of CRs is N CI and the typical region 
over which they have diffused is i?diff than the average CR density is n CI = 
37V cr /(47ri?j iff ). The pion decay luminosity scales approximately as 

L n o cx / n CT n R (R)4TrR 2 dR. (39) 
Jo 

For a uniform background density this means that 
L^o oc N CI n H , 

i.e. the total luminosity does not depend on i?diff ■ But for a supernova rem- 
nant evolving inside a stellar wind we have 

L^o cx 3N CI M w /{v w Rl iB ). 

Since i?diff oc Vi, this means that for core-collapse supernova remnants, the 
older the supernova remnant, the weaker the pion decay luminosity. This 
simple model does, of course, assume that the highest energy CRs are ac- 
celerated rather early in the life of a supernova remnant. For a supernova 
remnant inside a stellar wind that may well be the case, because the density 
and velocity are very high early on, and all the flux of particles entering the 



shock is higher in the beginning (see the case of SN 1993J in Sect. 2.3). This 
is in contrast to the case of a uniform medium, where the flux of particles 
entering the shock is slowly increasing with timej^] This difference in accel- 
eration properties for supernova remnants in different media will likely lead 
to larger maximum proton energies in core collapse supernova remnants, es- 
pecially those expanding in dense winds, and Type la supernova remnants 



Ptuskin and Zirakashvili 2005 Schure et al. 2010 \ 



If escape is indeed important, even for young supernova remnants, this im- 
plies that perhaps for Tycho's supernova remnant some of the 7-ray emission 
may not come from the shell, but from the vicinity of the supernova remnant. 
Indeed, there is a hint of an offset between Tycho's supernova remnant and 



8 For example, for the Sedov solution in a uniform medium, R oc t 2 ^ 5 we have 
F — 4nnn(R)R 2 v oc t 1 / 5 . Whereas for the Sedov solution inside a stellar wind, 
R oc t 2/3 , we get F = 4TTn H (R)R 2 v oc v oc t~ 1/3 ■ 
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the TeV 7-ray emission in the direction of a molecular cloud (Acciari et al. 



2011). 



The above analysis also implies that supernova remnants in a uniform 
density medium should still be visible for a relatively long time, as the total 
pion-decay luminosity is more or less constant. But in practice, at some mo- 
ment the total emission may come from such a large region that the 7-ray 
surface brightness will become low, and the supernova remnant CR "halo" 
will be difficult to detect, due to lack of contrast with the background emis- 
sion. Moreover, at large distance from the supernova remnant the diffusion 
coefficient is likely to become larger, as the CR generated, magnetic turbu- 
lence will be reduced if the CRs are spread out. Furthermore, the circum- 
stellar media of supernovae are usually more complicated than the simple 
models presented here. In particular, a stellar wind zone finishes in a region 
of enhanced density, the wind shell. This shell, once reached by escaping CRs. 
may be a source of enhanced pion decay emission. 

Over the next five years we will likely increase our understanding of 7-ray 
emission of young supernova remnants, as the current 7-ray observatories 
will continue their exploration of the Galaxy. But then the next big step 



forward will likely be made with the Cherenkov Telescope Array (CTA: CTA 



Consortium 2010 1, which will be ten times more sensitive, and have a broader 



energy range than current Cherenkov telescopes. This may provide more 
information on the importance of escaping CRs, and how they diffuse into 
and interact with the circumstellar medium, before becoming part of the 
overall population of Galactic CRs. 



2.14 GeV and TeV 7-ray observations of mature SNRs 



The first evidence that old supernova remnants can be GeV emitters too, was 
provided by Esposito et al. ( 1996 ), based on EGRET data. However, the point 
spread function of EGRET did not allow for an unambiguous identification 
of the 7-ray source with supernova remnants. Currently, Fermi has been able 
to confirm the EGRET results and has added many other mature supernova 
remnants to the list of the 7-ray sources (Table 2.14). 

Their 7-ray spectra are generally characterized by a broken power law 
with the peak of the emission in the GeV regime (e.g. Li and Chen 2011). 
The break photon energy is typically a few GeV and spectra above the break 
are very steep in comparison to the spectra as predicted by diffusive shock 
acceleration theory. For this reason, the GeV 7-ray telescopes as Fermi and 
AGILE are essential for detecting these intermediate aged remnants. Re- 
markably, shock waves of almost all these old remnants are interacting with 
dense molecular clouds, as evidenced by detection of hydroxyl (OH) masers 



at 1720 MHz (Wardle and Yusef-Zadeh 2002 Hewitt et al. 20091. The origin 



of the GeV 7-rays can be interpreted as hadronic processes (e.g 



Abdo et al. 



2010a). The high density in and near the supernova remnant suggests that 



the emission is likely due to pion decay, but bremsstrahlung cannot in all 
cases be ruled out, as this also scales with «h (e.g., Abdo et al. 2010d and 
section 2.13). However, new A GILE observations of two supernova remnants 
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Name 


Other 


GeV 


TeV 


Interacting? 




name 


7-rays 


7-rays 




G6.4-0.1 


W28 


Y 1 


Y« 


OH maser 3 


G8.7-0.1 


W30 


Y 4 


N 


OH maser 3 


G15.4+0.1 




N 


Y 5 




G22.7-0.2 




Y 8 


N 


N 


G23.3-0.3 


W41 


Y 7 


Y 9 


CO 8 


G31.9+0.0 


3C 391 


Y 4 


N 


OH-maser 3 


G34.7-0.4 


W44 


Y 9 


N 


OH-maser 10 


G43.3-0.2 


W49B 


Y n 


yl2 


H2 emission 13 


G49.2-0.7 


W51C 


Y 14 


ylS 


OH-maser 3 


G74.0-8.5 


Cygnus Loop 


yl5 


N 


CO emission 16 


G180.0-1.7 


S147 


Y 17 


N 




G189.1+3.0 


IC 443 


Y 18 


Y 19 


OH-maser 20 


G304.6+0.1 


Kes 17 


Y 21 


N 


OH-maser 22 


G318.2+0.1 




N 


maybe 23 


12 CO emission 


G348.5+0.1 


CTB 37A 


Y 4 


Y 24 


OH-maser 22 


G349.7+0.2 




Y 4 


N 


12 CO emission 


G359.1-0.5 


HESS J1745-303 


N 


Y 2 6 


OH-maser 27 



12 :i 



2r, 



Table 4 Table with mat ure supernova remnants that have been detected at GeV 



and/or TeV 7 energies 

fl2008c|), ^Hewitt et al 



et al. 



Abdo ct al. 
( 20091 



l2010a|);|Giuliani et alT] (|2010|), 2 |Aharonian 



Castro and Slanc (2010) 7jHofverberg et al 

(|2011bt °|Lafron et al.| <|2oTlb, f |Mehault et al.| (|2010|), B |Tian et al.| (|2007|), lAbdo 
et al.| (2010d J); |Giulian i et a l.| |2bll|l, lu |Hoffman et al.| J2005|), l1 |Abdo et al 



(|2010cD, 1Z |Brun et al.l ( |201lT 
Fiasson et al.| fl2008|), J& |Katagiri et al.|(|2011 



Keohane et al.j (12007 

T7 



ran 

>7P 



Abdo et al 



(2009) 



Katsuta et al. (2012) 



IN 



2011a| , ^ |Aharonian et al.| ^OOSaT^ 5 |Dubner et al. 



Abdo 



et al.|(|2010eT ; ITavani et al.| J2010| ) la |Albert et al.|(|2007a|);|Acciari et al.j J2009|) 



Hewitt et aL|(|2006|), z jWu et al.| J2011b , ^|Frail et al.l ( |1996|) , ^ |Hofverberg et aL 



2008b I, Wardlc and Yuscf-Zadch (2002 1. For a short descripion of most of these 
sources, see |Li and Uhen l 2Ull[ ). 



Aharonian et al 



(W44, W28 Giuliani et al. 2011 2010) reveal the characteristic low energy 



cut-off expected for pion decay (Figure [13| 

Several models exp lain the break and the s teep s pectra in terms of a 
crushed cloud model (|Aharonian and Atoyan 



1996 



Malkov et al. |2011 Uchiyama et al. 2010 Inoue et al 



Ohira et al. 2011a 



2010p . In this model, 



the GeV and TeV 7-ray emission is explained by pion decay from accelerated 
protons after the supernova remnant interacts with a high density molecular 
cloud. 

As CRs escape from remnants, they have a spectrum with a power-law 
index of appro ximately two (Ptuskin and Zirakashvili 2005). Aharonian and 
Atoyan| |l996) showed that 7-ray spectra CRs after escaping the remnant 
get steeper because of energy-dependent diffusion. The diffusion length after 
escaping from the supernova remnant, i?diff, is given by 



R diS (E) = ^/^D lSM (E)At 



(40) 



where Z?ism (E) and At are the diffusion coefficient of ISM and the elapsed 
time, respectively. Let Qsnr oc E~ s and Dism(E) cx E s be the CR spectrum 
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Fig. 13 Broadband spectrum of supernova remnant W44 with the best-fit model 
(red line). Overplotted is a hadronic model, revealing the characteristic low cut-off 
expected for pion decay (yellow curve). The green curves show the contribution of 
sever al emission mechanisms from energetic electrons (Figure from |Giuliani et al.| 



2011 reproduced by permission of the AAS. 



in the supernova remnant and the diffusion coefficient of ISM, respectively. 
Then, the CR spectrum at a certain point within the diffusion length, /diff, 
is given by 



,M ' R 3 diS (E) 



(41) 



The factor of -R^iff can be interpreted as the dilution due to the diffusion. 
Higher energy CRs have a longer diffusion lengths and are more diluted, 
so that the spectrum of escaping CRs becomes steeper than that of the 
source. The interaction between these escaping CRs and molecular clouds 
can produce 7-rays outside supernova remnants and can explain the uniden- 
tified very- h igh-energy 7-ray sources ( Aharonian and Atoyan| 1996 Gabici 



et al. 2009). Interestingly, HESS and MAGIC have now revealed several 



TeV sources, near, but not coinciding with, mature supernova remnants. 



This appears to be the case for W28 (Aharonian et al. 2008c I, IC 443 (Al- 



bert et al. 2007a) and G35. 6-0.4 (Torres et al. 2011). The idea is that the 



sources are molecular clouds that are illuminated as CRs that escaped from 
the supernova remnant are interacting with the dense material in the molec- 
ular cloud. Escaping CR electrons also produce 7-rays outside supernova 
remnants (Ohira et al. 2011b). 



In summary, it is unclear whether the GeV and TeV 7-ray emission from 
the youngest remnants is of hadronic of leptonic origin. In contrast, the GeV 
(and TeV) 7-ray emission from middle-aged supernova remnants is most 
likely explained by the hadronic scenario, giving credence to the paradigm 
that supernova remnants are capable of accelerating protons as well as elec- 
trons. As these old remnants are expected to no longer accelerate particles 
up to the knee, this does not provide evidence for supernova remnants being 
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able to accelerate particles to the knee. But since the new GeV/TeV obser- 
vations provide clear evidence for escaping CRs, acceleration up to or even 
beyond the knee is certainly possible. 



3 Balmer dominated Shocks 



A fast non-radiative shock in ionized gas produces none of the usual nebular 
emission lines, but if the gas is partly neutral, then some H atoms will be 
excited before they are ionized. They produce faint Balmer line emission fila- 
ments that trace the edges of several young Type I supernova remnants, along 
with a few older remnants and pulsar wind bow shocks. The Ha line profiles 



show broad and narrow components (Kirshner and Chevalier 1978), which 



provide valuable diagnostics for the proton and electron temperatures in the 
shocked gas, alon g with a means of detecting a shock precursor ( [Raymond] 



1991 Heng 2010). 



The neutrals pass through the shock front completely unaffected by the 
magnetic and electric fields in the shock or by plasma turbulence. They find 
themselves immersed in hot plasma moving at V s (l — 1/r), where V s is the 
shock speed and r is the compression ratio. One of three things can happen: 
they can be ionized by collisions with electrons or ions; they can be excited 
by electrons or ions, or they can undergo charge transfer with protons. Even- 
tually they all become ionized. All the H I emission is produced before the 
hydrogen is ionized, so it reflects conditions immediately behind the shock. 

Particles that are excited before any other interaction occurs produce 
Balmer line emission with the same velocity distribution that they had up- 
stream, and this is called the Narrow Component of the profile. Since the 
upstream plasma must be cooler than about 10,000 K in order to have a sub- 
stantial neutral fraction, one expects a line width below 25 km s _1 (FWHM). 
Other particles undergo charge transfer. This produces a population of neu- 
trals with a velocity distribution similar to that of the post-shock protons, 
but weighted by the velocity times the charge transfer cross section. In shocks 
faster than about 1000 km s _1 , the charge transfer can go into an excited 
state and produce an Ha photon. Either charge transfer to the excited states 
or collisional excitation of the fast neutrals will produce a Broad Component, 
whose line width is on the order of the post-shock proton thermal speed 



(Chevalier and Raymond 


1978 


Chevalier et al. 


1980 


2007 


Heng et al. 


2007 van Ade 


sberg et al. 


2008 


!• 



Heng and McCray 



The line widths of the Narrow and Broad components give the pre- and 
post-shock temperatures fairly directly. The intensity ratio of the Broad and 
Narrow components depends on the ratio of ionization and charge transfer 
rates, so it provides a means of determining the electron-to-ion temperature 
ratio immediately behind the shock. In the following, we discuss the impli- 
cations for shock precursors and for electron-ion temperature equilibration. 
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Fig. 14 Ha line emission of the southwest shock of supernova remnant 0509-67.5 in 
the Large Magellanic Cloud. The solid line indicates a fit to the line, with a broad 
and a narrow gaussian line profile convolved with the spectral resolution. Note that 
the width of the narrow component is determined by the spectral resolution. Figure 
was published in Helder et al.| ( |2010| , reproduced by permission of the AAS. 



3.1 Anisotropy and non-Maxwellian distributions 



Nearly all analyses of X-ray, optical and UV emission from shocks assume 
isotropic, Maxwellian velocity distributions, but in a collisionless plasma that 
may not be the case. Several types of complex distributions are possible. 



Neutrals following charge transfer: Heng and McCray (2007), Heng et al. 
( |2007[ ), |van Adelsberg et al.| ( |2008[ ) and |de Laat| ( |2011| ) have carried out 
detailed calculations of the H I velocity distribution in the shocked plasma. 
For moderate shock speeds, the distribution is nearly Maxwellian, but at 
speeds above about 2000 km s _1 , the rapid dropoff of the charge transfer 
cross section with speed produces distributions that are narrower in the flow 
direction than the perpendicular direction, and they are not Maxwellian. 



Pickup Ions: Ohira et al. (2009b) and Raymond et al. (2008) discuss the 



velocity distribution of protons produced when H atoms that pass through 
the shock as neutrals become ionized downstream. Because they suddenly 
find themselves immersed in a rapidily moving magnetized plasma, they be- 



have likepickup ions in the solar wind (Moebius et al. 1985 Gloeckler et al. 



1993 Isenberg 1995). They initially have the form of a monoenergetic ring 



beam gyrating around the magnetic field. That unstable distribution gener- 
ates plasma waves and relaxes to a bispherical distribution symmetric about 
the magnetic field direction. Ohira et al. ( |2009b[ ) discuss the wave genera- 
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tion, magnetic field amplification and shock structure modification that can 



result, and those changes could affect CR acceleration. Raymond et al. (2008) 
predict the signatures that could be detectable in the Ha line profile. 

Power-law tails: cosmic-ray acceleration implies a power-law tail at some 
level, and velocity distributions in the solar wind are often seen to have a 
K-distribution character, consisting of a nearly Maxwellian core with a power- 
law tail. The power-law tail must contain a large fraction of the total energy 
in order to produce detectable wings on the Ha profile. 

Observations of Ha line profiles are generally interpreted by fitting Gaus- 



sians for the Broad and Narrow components, but Raymond et al. (2010) 
show the profile from a shock in Tycho's supernova remnant that cannot be 
fit with two Gaussians. They are not able to discriminate among the pos- 
sibilities that the profile results from pickup ions, a k distribution, a hot 
precursor, or merely the superposition of different shocks within the spectro- 
graph slit. Analysis of other spectra from Tycho may help to sort out the 
options. 



3.2 Shock precursors 



The diffusive shock acceleration model for acceleration of CRs in shock waves 
requires a precursor region where energetic particles streaming away from the 
shock generate turbulence, and where the turbulence scatters par ticles back 
toward the shock s o that first order Fermi acceleration can occur ( Blandford 
and Eichler 19871) . The plasma in the precursor is compressed, heated and 



accelerated toward the shock speed (Boulares and Cox 1988 Vladimirov 



et al. 20081. The scale length of the precursor is given by the diffusion 



coefficient of the most energetic cosmic rays, D (see equation 12), divided by 
V s , and it should be on the order of an arcsecond for typical young Galactic 
supernova remnants. A precursor unrelated to CRs is also possible, in that a 
few percent of the neutrals in the Broad Component can overtake the shock, 
and if they deposit energy and momentum there they can dece lerate (in the 
shock frame) and heat the incoming gas ( Hester et al.| 1994 Smith et al. 
1991). Its length scale is the charge transfer length scale, which is typically 



0.1 to 1 arcsecond. 

Indirect evidence for a shock precursor comes from the narrow compo- 



nent widths measured in 9 supernova remnants ( Sollerman et al. 2003 ). With 
the exception of SN1006, the widths are 35 to 58 km s _1 . At the tempera- 
tures implied by those widths, hydrogen would be fully ionized, so the heat- 
ing must occur in a precursor narrow enough that much of the neutral gas 
passes through the precursor before being ionized. This implies thicknesses 
compatible with th e CR precursor, which in turn would imply D values o f 



10 25 cm 2 /s (Hester et al 



1994 



1991 



2007). 



10 2 

The one supernova remnant that shows no enhanced line width is SN 1006 



Smith et al. 



Lee et al. 



(Sollerman et al. 2003]). This could be related to the lack of non-thermal 



emission in that section of the supernova remnant, but it might also be that 
the density is so low that neutrals pass through the precursor without inter- 
acting at all, so that their line width does not reflect heating in the precursor. 
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More direct evidence for a precursor comes from the observation of faint 
emission ahead of the main shock filaments in Tycho's supernova remnant. 



Spectra from Subaru show this to be narrow component emission (Lee et al. 



2007), and an HST image shows that it trails off gradually ahead of the shock 
(|Lee et aT| 20101. By combining the observations with models by |Wagnerj 



et al. 



Lee et al. (2010) obtained a precursor length scale of about 



2009f;_ 

D cm and a peak temperature in the precursor of 80,000-100,000 K. 

1 The neutrals may affect D, since they could 



6 x 10 

This implies D 

damp the waves that scatter energetic particles ( |Drury et al. 1996), which 
lowers the maximum energy the CRs can achieve and reduces the overall 
acceleration efficiency. 

The interaction between neutrals and a CR precursor may be much more 
complex than is suggested by the estimates of scale length and temperature. 
A strong precursor accelerates the pre-shock plasma to a substantial fraction 
of the post-shock speed. When a neutral becomes ionized in the precursor, 
either through charge transfer or collisional ionization, it behaves much like 
a pickup ion in the solar wind (|Ohira and Takahara 2010 Raymond et al 



2011). The resulting population of protons may retain an anisotropic, high 



speed velocity distribution, or it may come into equilibrium with the ther- 
mal protons by way of plasma wave interactions. In either case, if the neutral 
fraction is substantial there can be important effects on the heating and com- 
pression in the precursor and on the injection of particles into the acceleration 
process. 

It should be kept in mind that Broad component neutrals that overtake 
the shock, as suggested by Hester et al. ( 1994 1 and Smith et al. ( 1991 ) might 
also account for the observed Narrow comp onent width of Ha and for the Ha 
ahead of the main shock front observed by Lee et al. ( 2007[ 2010). Complete 
models are not yet available, though Lim and Raga (1996) found that this 
process had little effect on the Ha profile. However, they did not consider 
the possibility that the protons formed from fast neutrals that penetrate 
the shock deposit a significant fraction of their energy by generating plasma 
waves. 



3.3 Temperature versus shock speed 



Once a shock deposits a substantial fraction of its kinetic energy into acceler- 
ating particles, the structure of the shock will change. A particle accelerating 
shock will have a lower post-shock temperature that a non-accelerating shock. 
In this section, we first describe how the post-shock temperature relates to 
the post-shock cosmic-ray pressure and cosmic-ray energy flux escaping the 
system (3.3.1). In 3.3.2 we describe how the proton and electron tempera- 



tures can differ and where this causes caveats in determining the cosmic-ray 
acceleration efficiency. Section [3.3.3| lists the observations that have been car- 
ried out to determine cosmic-ray acceleration efficiency parameters based on 
post-shock plasma temperatures. Section |3.3.4| concludes with the relation 
between the post-shock pressure and the cosmic-ray energy flux escaping the 
system. 
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3. 3. 1 Physical parameters 

If a shock accelerates CRs, it changes the physical parameters of the plasma 
behind the shock front. Firstly, the CR pressure increases as the shock is ac- 
celerating more efficiently. The total pressure behind the shock front remains 
approximately constant, implying that the local thermal pressure, and there- 
with the temperature, drops. Secondly, the effective equation-of-state (7 S ) of 
the post-shock plasma changes from 7 S = 5/3 (non-relativistic) t o 7 S = 4/3 
(relat ivistic) if the post-shock pressure gets more CR dominated ( Chevalier] 



1983) 



7s = — = (42) 
7 S - 1 3(1 + w C r) 

Here, u>cr is the fraction of the post-shock pressure that is contributed by 
CRs. Using the conservation of mass, momentum and energy over the shock 
front, and taking into account that a fraction (ccr) of the energy flux through 
the shock {\pVg) is escaping the system, we get the following expression for 



the compression ratio over the shock front (c.f., Bcrczhko and Ellison 1999 



Vink 2008a): 



_ G +y/(P-{l-e CK ){2G-l) 
r l- ecR ' 

where G = §ii>cr + §■ Altogether, this leads to a mean post-shock tem- 
perature (T) of 



(kT) = (1 - Wc r)-(1 - l/r)^m p V s 2 . (44) 

Here, /i is the mean particle mass in units of the proton mass (rn p ), which 
is ^0.6 for solar abundances. This has to be compared to the temperature 
behind non-accelerating shocks: (kT) — 3 /i6fj,m p Vg. To quantify the effect of 
CR acceleration, we define a parameter 

/3 = ( feT >/3/16 M mpK 2 - 

So, for a non-accelerating shock, f3 — 1 and a lower /? implies more efficient 
shock accelerating CRs. 

3.3.2 Electron-proton temperature equilibration 

To measure the cosmic-ray acceleration efficiency from the temperature deficit, 
we have to know the mean plasma temperature. Most young supernova rem- 
nants emit thermal X-ray emission, enabling us to determine their electron 
temperatures from the shape of the bremsstrahlung continuum. Unfortu- 
nately, directly behind the shock front, the electron temperature is not nec- 
essarily equal to the temperature of other atomic species. They might even 
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Fig. 15 /3 (black line) and compression ratio (grey dashed line) over 
front as function of energy flux escaping the sho ck (as fract i on of the 
energy, Jtot = |pK 3 )- Based on the equations of Vink et al. (20101. 



the shock 
incoming 



differ as much as their mass ratio (a factor of 1836 for electrons to protons). 
Therefore, the electron temperature at the shock front is a poor indicator for 
the mean plasma temperature. 

In time, the temperatures will evolve to the mean plasma temperature. 
If Coulomb interactions are the only process contribution to equilibration, 
this will happen after n c t ~ 10 12 s cm~ 3 , where n c is the electron den- 
sity, and s the time in seconds (Vink 2004). However, it is rather uncer- 



tain at which electron to proton temperature ratio a plasma has if it was 
just shocked. From optical observations, a relation of was proposed for 
V s > 400 km s -1 , below 400 km s -1 the electron and proton temperature 



were similar (Ghavamian et al. 2007[ ) (however, see van Adelsberg et al. 



2008 



Helder et al. 



2011 1. This change of behavior at 400 km s 1 might be 



explained if one realizes that this if very similar to the electron velocities in a 
neutral medium with a temperature of 10 4 K. Th erefore , the electrons might 
not feel the shock as a discontinuity (Bykov and Uvarov, 1999 Bykov, 2004 1. 

It is not well-understood how cosmic-ray acceleration affects ion-electron 
temperature equilibration. One possibility is that the heating of electrons in 
the CR precursor in combination with the lower Mach number will lead to 
temperature equilibration for shocks with velocities > 400 km s" 1 . Addition- 
ally, Patnaude et al. ( |2009 ) argue that the increased compression ratio of 
particle accelerating shocks (section [4 ) accelerates the equilibration, as the 
ionization time (n e t) is a function of the density too. 



3.3.3 Observations 



Hughes et al. (2000) measured the electron temperature behind the shock 
front of SNR 0102.2-7210 in the small Magellanic Cloud as well as the shock 
velocity based on X-ray data. They found that, even after correcting for 
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temperature equilibration effects, the measured electron temperature was 
still too low for the measured shock velocity, indicating that SNR 0102.2- 
7210 has efficiently acce lerating shocks. 

Salvesen et al. ( 2009 ) did a similar study for the ~ 10,000 year old Cygnus 



Loop remnant, based on electron temperatures obtained from X-ray spectra 
as extracted from ROSAT-PSPC data. The shock velocity was determined 
from proper motion measurements taken from optical images obtained with 
the Palomar Observatory Sky Survey in two epochs, 39.1 years apart. From 
this study follows that the post-shock cosmic-ray pressure is negligible. 



Helder et al. ( 2009 ) did a similar study on the TeV 7-ray remnant RCW 86 
( Aharonian et al. |2009[ ) at the location where the shock front shows X- 
ray synchrotron emission. They used the Ha line to determine the post- 
shock proton temperature, and 2 epochs of Chandra data to measure the 
proper motion of the shock front. The shock velocity was significantly higher 
than one would expect from the post-shock proton temperature, resulting in 
(3 = 0.3. 

However, one has to caution that for the latter two studies, the proper 
motions were determined using emission from a different wavelength than 
the emission used to determine post-shock plasma temperature. This has the 
caveat that one might match physical parameters for shocks that happen to 
be in the same line of sight, but are not part of the same physical system. 



Helder et al. (20101 measured the post-shock proton temperature as ob- 
tained with VLT/FORS2 for the LMC remnant SNR 0509-67.5 (Fig. fill. 



They combined this with a shock velocity determined from high resolu- 
tion X-ray spec t ra as obtained with the XMM-Newton RGS instrument 



(Kosenko et al. 2008). The resulting shock velocities were 5000 km s 1 



and 6000 km s _1 for the SW and NE shocks respectively. These velocities 
seem to be confirmed by results from a preliminary study carried out by 
Hovey and Hughes (2012), as they measure an average shock velocity of 6500 
±200 km s _i around the remnant, based on a proper motion study using 
two ep ochs of Hubble Space Telescope data. The velocities of |Helder et al. 
(120101) resulted in a value of f) = 0.7 for the SW shock. The NE shock does 



not show evidence for cosmic-ray acceleration, but note that the uncertainty 
on the post-shock proton temperature is rather large, and therefore efficient 
cosmic-ray acceleration is not excluded either. 



3.3.4 Relation wcr and ecR 



Equation 44 shows that measuring the fraction by which the temperature 



is lowered f3 in itself docs not provide information on the conditions of the 
cosmic-ray pressure of cosmic-ray escape of the shock. To solve this, |Helder| 



et al. (20091 and 



leading to lower 



imits for wcr'- 



Helder et al. (2010) used a curve as an upper limit for ecR, 



£cr 
w C R 



(1 - l/rf 



(45) 



This curve is based on the asymptotic particle spectra of a strong shock 
Malkov (1999). In a more recent study, Vink et al. (2010) calculated the 
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relation between e^R, and wcr, using the standard Rankine-Hugoniot jump 
conditions over a shock, this time including a cosmic-ray precursor. Surpris- 
ingly, the ecR, — wcR-relation resulting from this study overlaps exactly with 
the upper limit in equation 45 for a 7 of the CRs of 5/3. For a 7 of the CRs 

0.5 and e CR = 0.2 
0.03. 



of 4/3, we get the following results: for RCW 86, Wcr 
and for the SW shock of SNR 0509-67.5, w C r = 0.21 and e C R 



4 Increased compression ratio 



As soon as a shock becomes efficient in accelerating, the density behind the 
shock front will increase for two reasons. The equation of state of the post- 
shock plasma will soften, leading to a higher overall compression ratio. Also, if 
a shock loses its energy in the form of CRs escaping, it compensates for these 
energy losses by making the plasma behind the shock front denser, in order to 
maintain pressure equilibrium (equation 43 and Fig. 15). This immediately 



implies that the forward shock will be closer to the reverse shock and contact 



discontinuity than for a non-accelerating shock ( Decourchelle et al. 2000 



Ellison et al. 


2004 


Kosenko et al. 


2011) 



2011 ). The grey line in Fig. 15 shows the 



relation of the compression ratio as function of ccr- 

This provides a way for measuring such an increased density. If one can 
accurately trace the outer edges of the shocked ejecta (i.e., the contact discon- 
tinuity), this can be compared with standard, non-CR accelerating models 
of supernova remnant evolution. This was done for the first time by | Warren 
eFaLl ( p005| for the Tycho supernova remnant using a principal component 
analysis to determine the location of the forward shock, contact discontinu- 
ity and reverse shock. Indeed, in Tycho, the contact discontinuity appears to 
be closer to the shock front than expected from standard models, implying 
efficient cosmic-ray acceleration ( |Kosenko et al. 2011| . 

In interpreting the outer edges of the ejecta in terms of the contact dis- 
continuity, one should be careful. Hydrodynamical instabilities can bring the 



eject a in fron t of the contact discon tinuity (Blondin and Ellison 2001 Wang 
2011). Cassam-Chena'i et al. (2008) did a similar study on the SN1006 rem- 
nant, and revealed that the ejecta are closer to the forward shock at locations 
coinciding with X-ray synchrotron emission and TeV 7-ray emission. How- 
ever, they found the contact discontinuity to be too close to the shock front, 
for efficiently accelerating models, even for the regions that did not coin- 
cide with TeV 7-ray emission of X-ray synchrotron emission. In a subsequent 
study, Miceli et al. (20091 found that the ejecta are less close to the shock 



front. Also, their study showed the ejecta to be similarly close to the shock 
front at the locations of the TeV 7-ray emission as well as at other shocks. 
In fact, they concluded the shocks to be modified by cosmic-ray acceleration 
everywhere. 

Note that the separation between shocked ejecta and the forward shocks 
might also be reduced if the forward shock encounters a density jump in the 
ambient medium (see e.g. Kosenko et al. 20101. 
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4.1 Effect of CRs on Supernova Ejecta Morphology and Clumping 



Blondin and Ellison (2001) carried out hydrodynamical simulations on the 



evolution of supernova remnants, including the effects of particle acceleration. 
In particular, they studied the effect of high compression ratios resulting 
from efficient cosmic-ray acceleration. They found a decrease in the distance 
between the forward and contact discontinuity with respect to a shock that 
was not accelerating particles. Moreover, in cases where particle acceleration 
is particularly efficient, instabilities at the contact discontinuity might reach 
the forward shock . This will cause the shock front to be perturbed by these 
instabilities. Rakowski et al. (2011 ) found protrusions along the eastern rim of 



SN 1006. This was interpreted in the context of an upstream magnetic field, 



amplified by streaming CRs Bell (2004 
instabilities post shock. 



2005), enhancing Rayleigh- Taylor 



4.2 Fast Moving Supernova Ejecta Fragments 



Aschenbach et al. ( 1995 ) discovered protruding structures in the Vela su- 



pernova remnant, using the ROSAT satellite and dubbed them 'shrapnels'. 
Miyata et al. (2001 ) used Chandra to investigate 'shrapnel A' and found that 



its silicon abundance is above solar. Cassiopeia A also has fragments of ejecta 



outside of the main rim. These were first detected in optical (e.g. van den 
Bergh 1971 Kirshner and Chevalier 1977 and the references therein). These 



fragments were interpreted in the context of shock- wave emission models by 
Raymond (1976), resulting in the conclusion that their abundances are con- 
sistent with uncon taminated material from the core of a massive star. [Laming] 
and Hwang ( 2003 ) analyzed several knots of Cassiopeia A as observed with 
Chandra. 



5 Molecular ion diagnostics of CR acceleration sources 



For a long time, the ionization of the interstellar medium by CRs has been 
considered an important factor in the phases of the interstellar medium ( Field 



et al. 1969 Spitzer 1978 1. The chemical reactions that take place in molecu- 



lar clouds are determined by CR ionizat i on rat es (e.g. Herbst and Klemperer 
Dalgarno] 2006 ) . |Padovani et al.| ( 2009 ) calculated the gas ionization 
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rates in molecular clouds assuming a CR spectrum extrapolated to low ener- 
gies from observed CR proton and electron spectra. Hereby, they took into 
account CR particle energy loss processes. Their model agreed with obser- 
vational data, only if the energy flux of the CR spectrum increases at low 
energies (below ^100 McV). The diffusion of CRs within molecular clouds are 
thought to be the dominated by the magnetic field orientation and strength. 
Padovani and Galli (2011) modeled this diffusion and found that magnetic 



mirroring is more important than magnetic focussing, resulting in a reduc- 
tion of the CR ionization rate with a factor of about 2-3 in the core of the 
molecular cloud. 
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Molecular ion line observations may serve as a diagnostics of processes in 
CR acceleration sources. Near the IC443 supernova re mnant, 



a column 



density of 3x 10 cm s was found in two sight lines ( Indriolo et al. , 2010 1 . 



From these measurement, the ionization rate was found to exceed that of the 
average ionization rate of molecular clouds. However, the other four sight 
lines were consistent with typical Galactic values. A similar study on the 
W51C supernova remnant was carried out by Ceccarelli et al. (2011). Also 



for this remnant, the cloud ionization degree is highly enhanced, resulting 



T0" 



-15 



s . This is about 100 times the average 



in a CR ionization rate of 
value for molecular clouds. 

Modeling of the ionization rate in interstellar material around a CR accel- 
erating source shows that in infrared, the strongest lines of H^" are between 
3.0 and 4.6/xm and H j~ line s of similar intensity are visible between 4.0 and 
6.0 /xm |Becker et al.| ( |2011| ). 

The Herschel key proj ect "Chemica l HErschel Surveys of Star forming 
region" (CHESS) (see e.g. Ceccarelli et al. 2010) collected high resolution 



molecular line spectra of star-forming regions that can be used to constrain 
the gas ionization rates in the molecular clouds. Also, Gerin et al. (2010) 



reported the detection of the ground state doublet of the methylidyne radical 
CH at 532 GHz and 536 GHz with the Herschel/HIFI instrument along the 
sight-line to the massive star-forming regions G10.6-0.4 (W31C), W49N, and 
W51 where some CR acceleration sources are residing. 



6 Particle acceleration in the Galactic Central regions 



The central 200 pc of the Galaxy are intense high energy emitters, from 
X-rays to TeV 7-rays. The physical conditions in these regions are indeed 
exceptional and make the Galactic Center (GC) a special place as regards 
particle acceleration. Bright non-thermal radio filaments trace magnetic fields 
of up to a mG ( Ferriere 2009 ) an d large scale non-the rmal emission suggest 
the average field is at least 50 /xG (Crocker et al. 2010[ ). The Cent ral Molecu- 
lar Zone conta ins about 50 million solar masses of molecular gas ( |Morris and 
Serabyn 1996 ) in which physical conditions are quite unusual, with high tem- 



peratures and a large velocity dispersion. The abundance and high density 
of the gas are such that intense star formation is taking place in this re- 
gion. Three of the most massive young star clusters of the Galaxy lie within 
30 pc of the GC: the Arches, the Quintuplet and the central cluster. The 
Quintuplet harbors what could be the most massive star in the Galaxy: the 
Pistol star, a Luminous Blue Variable candidate with an estimated initial 
mass of 200 M (Yungelson et al. 2008[). The Arches cluster is one of the 



most massive and dense star clusters in the Galaxy. Its age is just 1-2 My r 



and its central density reaches 10 5 M© pc d (Portegies Zwart et al. 2010). 



The inner 10 pc are dominated by a bright supernova remnant, called Sgr A 
East. The radio shell is filled with unusually hot thermal X-ray emission and 
is consistent with a single star explosion about 10000 years ago. Finally at 
the very center, the faint source Sgr A*, is powered by a supermassive black 



hole of 4 10 M (Genzel et al. 2010 and references therein). The black 



hole is particularly quiet nowadays with a quiescent bolometric luminosity 
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of about 10 36 erg/s but there is good evidence that it has not always been 
so dim. 

This exceptional environment and these extreme conditions make the GC 
region an interesting environment for testing particle acceleration models. 
After this very brief description of the scenery in the central 200 pc, we review 
some observational evidence for particle acceleration in the Galactic Central 
regions and discuss their possible origins. We discuss the observation of non- 
thermal emission from giant molecular clouds in the GC, first in non-thermal 
hard X-rays, possibly produced by sub-relativistic particles, and second in 
the TeV range, tracing an excess of very high energy CRs. Finally, we discuss 
the recent discovery of large bubbles of GeV emission extending above and 
below the GC. 



6.1 X-ray emission from GC molecular clouds 



One of the great legacies of the large modern X-ray observatories is the 
detailed mapping of the GC region. Among the main results, is the discovery 
that Sgr A* is both remarkably faint nowadays with a quiescent luminosity 



of 10 33 erg/s (|Baganoff et al 
(Porquet et al. 



20031 and subject to frequent X-ray flares 
2008[ ), whose peak luminosities reach only a few 10 35 erg/s 
or 10 _y LEdd for a 4 x 10 6 M Q black hole. The surveys of the central 100 
pc have also discovered that several massive molecular clouds in the GC 
display intense Fe Ka 6.4 keV line emission (|Koyama et al. 1996 |2007 l 



and non-thermal continuum (Sunyaev et al. 19931 IRevnivtsev et al. 12004 



Bclanger et al. 2006). The origin of this emission has been strongly debated. 



The propagation of energetic charged particles inside the molecular clouds 
has been invoked as a possible explanation of the observed X-rays from the 
region. 

Two distinct processes can produce hard X-ray continuum in combination 
with neutral iron line emission. First, low-energy CR electrons, with energies 
up to a few hundred keV, diffusing into dense neutral matter will produce 



nonthermal bremsstrahlung and collisional K-shell ionization (Valinia et al 
2000). Second, subrelativistic ions propagating in the molecular cloud can 



radiate through inverse bremsstrahlung in the hard X-ray domain and create 

It is likely that the density of low 



2009) 



iron Ka vacancies (Dogiel et al. 
energy CRs is larger in the central regions than in the rest of the Galaxy. 
The presence of strong non-thermal radio filaments in the vicinity of several 
molecular clouds makes it reasonable to consider the presence of a large 
density of CR electrons there. This is also supported, as noted by |Yusef-Zadeh| 
et al. ( |2007 ), by the observation of d iffuse low-frequency radio emission in the 
central regions ( LaRosa et al. 2005 ) and the high estimates of the ionization 



rate compared to the values obtained in the Galactic disk (Oka et al. 2005 ). 



A possible mechanism to explain the production of such subrelativistic 
particles a nd t he associated emission has been proposed by (Bykov 2003 



see section 4.2) in terms of fast moving knots resulting from supernova ex- 
plosions. 



Bykov ( 2002 ) has modeled this phenomenon taking into account injection, 



transport, first and second order Fermi acceleration as well as Coulomb losses 
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which are dominant for low energy particles. The efficiency of bremsstrahlung 
radiation compared to ionization losses is more favorable in such enriched 
clumps which alleviates the energy requirement for such a process to take 
place. In particular, an O-dominated clump with an X -ray luminosi t y of 10 31 
erg/s will lose ~ 10 35 erg/s through Coulomb losses (Bykov et al. 2008a|. 
Unfortunately, the recent observations of variability of the 6.4 keV line 



emission in several distant clouds ( Muno et al. 2007 Inui et al. 2009 Terrier 



et al. 2010) is contradicting the idea that CRs are producing the 6.4 keV 



line emission. Besides, the observation of an apparent superluminal propaga- 
tion of the 6.4 keV emission in a molecular structure ( Ponti et al.[ 2010) is 
also a clear sign that the clouds are in fact X-Ray Reflection Nebulae. Their 
emission is the result of Compton scattering and K-shell photo-ionization of 
neutral or low-ionized iron atoms by an intense X-ray radiation field gener- 
ated in the GC. The required luminosity of illuminating the source is about 
10 39 erg/s making the supermassive black hole, Sgr A*, the most likely can- 
didate. 

Not all Fe Ka emission is expected to be variable though, as is visible 
for instance on most clouds surrounding the Arches cluster region whose 
line flux have stayed constant over more than 8 years of XMM-Newton ob- 



servations (Capelli et al. 2011). The best explanation for this (up-to-now) 



stationary emission is that it is produced by subrelativistic CR accelerated 
locally, possibly by collective effects in the cluster itself. 



6.2 An excess of very high energy CRs in the central 100 pc 



Besides the emission of the pulsar wind nebula GO. 9+0.1 and the GC point 
source, the Cerenkov telescope system HESS has observed diffuse emission 
over 200 GeV extending over the central degree. The spatial distribution of 
this emission is well correlated with that of interstellar matter which sug- 
gests that the very high energy photons are produced by CR ions interacting 
with dense clouds (Aharonian et al. 2006). The photon spectral index of 



this emission is 2.3, somewhat harder than the CR spectrum measured at 
Earth. This proves that the density of high energy CRs (E>1 TeV) in the 
central 100 pc is significantly larger than in a solar neighborhood. One or 
several CR accelerators must therefore be presently or very recently active 
there. Aharonian et al. (20061 have proposed that the supernova remnant 



Sgr A East is the source of these CRs. Taking into account the diffusion time 
of particles in the central regions, they showed that about 10% of the ki- 
netic energy of the supernova remnant could explain the observed emission. 
Time dependent escape of energetic particles from supernova remnants has 
been proposed and studied by Gabici et al. ( 2009 1 and probably observed in 



systems such as W28 ( Aharonian et al. 2008c ) 



The propagation of CRs in the complex environment of the GC is most 
likely quite complex. For instance, the large poloidal fields traced by non- 
thermal filaments might limit the diffusion of CRs in the Galactic plane. 
The idea of a single accelerator inj ecting CRs in the GC has therefore been 
questioned ( Wommer et al. 2008 see for instance) and distributed turbu- 



lent acceleration in the inter-cloud medium has been proposed by several 
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authors (Melia and Fatuzzo 2011 Amano et al. 2011). We note that the 



massive young star clusters in the region can also provide significant particle 
acceleration through the collective effects of stellar winds and shocks. The 
interaction of the supersonic winds of the stars leads to the creation of a 
large supcrbubble filled with a hot and tenuous plasma. The propagation of 
supernova shocks inside the bubble create regions of enhanced turbulence 



where effective particle acceleration can take place (Bykov 2001 Parizot 



|etlJ~[ [20041 |Higdon and Lingenfelter] [20051 |Ferrand and Marcowithj|2010[ , 
this aspect is further discussed in the next section). 

High-energy and very high energy 7-rays have been detected in regions 
spatially coincident with massive young clusters in the Galaxy: HESS has 



observed TeV photons possibly coming from the Westerlund 2 (Aharonian 
et al.[ |2007[ |HESS Collaboration et al.[ |2011b[ ) and Westerlund 1 clusters 



(HESS Collaboration et al. 201 lap . Very recently, an intense G eV emission 



pervading the Cygnus X region has also been found with Fermi (Ackermann 



et al. 2011 ). Although the nature and origin of these emission regions are not 



completely proven yet, they are arguing in favor of the collective effects inside 
massive star associations as effective particle accelerators. It is therefore very 
likely that the exceptional clusters at the GC are actively producing high 
energy particles, although future TeV observatories will be necessary to prove 
it. Thanks to their improved angular resolution compared to HESS they 
should be able to resolve the contribution of clusters to the diffuse GC TeV 
emission. We also note that the winds from the central cluster, containing two 
dozens of supergiant luminous and Wolf-Rayet stars as well as numerous O 
and B objects ( jPaumard et al. 2006), are most likely powering the accretion 
onto Sgr A* and might also be at the origin o f the strong point-like emi ssion 
observed by HESS at the position of the GC (Quataert and Loeb 2005). 



6.3 Large bubbles of 7-ray emission powered by the GC 



Recently, thanks to specific procedure of foreground subtraction with the 
Fermi data, Su et al. (2010) found a pair of large structures extending up 



to 50 degrees above and below the Galactic plane in the direction of the 
GC. These structures, called the "Fermi bubbles", are detected above a few 
GeV and have a hard spectrum extending up to 100 GeV. They see m to be 
associated to large scale structures observed in X-rays with ROSAT (Bland 



[Hawthorn and Cohen 2003[ ) and with an excess of radio emission found in 
WMAP data. At the distance of the GC, they have a luminosity of 4 x 10 37 
erg/s and extend up to 10 kpc outside of the Galactic dis k. The estimated 
energy content of the bubbles is of the order of 10 55 erg (Bland-Hawthorn 



and Cohen 2003 Su et al. 2010). 



The nature and the origin of the particles creating this 7-ray emission are 
still unclear. Inverse Compton emission of electrons requires only a fraction 
of the total energy in the form of energetic particles but requires in-situ 
acceleration to have very high energy electrons up to 10 kpc above the disk. 
Energetic protons or nuclei have much longer lifetimes and can be injected 
from the Galactic central regions but require a much more substantial fraction 
of the total energy in the bubbles. 



57 



The mechanism providing the large energy input is also quite uncertain. 



Su et al. (2010) discussed a possible starburst phase in the GC a few million 
years ago. Today, a more modest but constant injection of 10 39 erg/s of CRs 
from the GC to the halo is inferred from the 7-ray emission from the central 
molecular zone and its infrared luminosity (Crocker et al. 20111. If this 



injection is integrated over a billion years, it might provide the correct energy 



budget (Crocker and Aharonian 2011), although the particles have to be 



confined on extremely long timescales. Other mechanisms outline the possible 
role of the supermassive black hole. An intense AGN phase at Eddington 
luminosi ty accompanied by jet s or outflows a few millions years ago has been 



invoked ( jZubovas et al. 2011) as well as the r ecurrent (every 10 -10 5 years) 



accretion of stars captured by the black hole (Cheng et al. 2011 ) 



A lot remains to be understood on the physics of these large 7-ray bub- 
bles as well as on the processes producing the high energy emission in the 
molecular clouds of the central 200 pc. Nevertheless, they remind us that the 
GC plays a special role in the production of Galactic CRs. 



7 Alternatives for SNRs being main sources for GCRs 



7.1 Particle acceleration in superbubbles 

Somewhere between 10 17 eV and 10 19 eV, the CR origin is expected to transi 
tion from Galactic to extragalactic origin (e.g. Hillas 2005| |Aharonian et al. 
2011). For this reason, Galactic cosmic-ray sources probably accelerate cos- 



mic rays to above the knee. Important in models of CR acceleration is the en- 
vironment in which the supernova evolves. Core collapse SNe usually reside in 
associations of young massive stars and molecular clouds. These associations 
usually create superbubbles i.e., large caverns filled mainly with a mixture 
of hot and warm gas produced by multiple supernova explosions. These su- 
perbubbles can have a significant effect on the structure of the disk and halo 



of the Galaxy (e.g. Norman and Ikeuchi 1989 Heiles 1990). Norman and 



Ikeuchi ( |1989[ ) found that superbubbles during their outbursts create mass 



and energy flows from the disk into the halo and vice versa. The ongoing 
shocks marking the outer edges of a superbubbles are a plausible location of 
CR acceleration. Additionally, the material in the superbubbles possi bly has 
an enriched composition with Wolf-Rayet winds and SN ejecta (e.g. |Bykov 
and Fleishman] [1992} |Axford[ [19941 |Klepach et aTl[2^0^|Bykov||200^l(|Birms 



et al. 2007 Ferrand and Marcowith 2010 1 



Observational studies on OB-associations and young globular cluster have 
been carried out both in the Galaxy and the Large Magellanic Cloud. A 
detaile d study on the Galactic OB-association Cygnus OB-2 by |Knodlseder 
(120001) showed that it has 2600 ± 400 OB members, out of which 120 ± 



20 are O stars. The corresponding high number of X-ray sources (Wright 
and Drakel 120091) suggests that Cygnus OB2 is one of the most massive star 



forming regions in the Galaxy. Given its relative small size, it is inevitable 
that the region contains strong winds closely together, accumulating in strong 
shocks. These are the regions that are particularly expected to be particle 
accelerators (e.g. Eichler and Usov 1993 Torres et al. [12004"] ). 



58 



Indeed, superbubble 3 Dor C in the Larg e Magellanic Cloud appears to 
emit non-thermal X-rays ( Bamba et al. 2004 ) . The shape of this non-thermal 
emission, coincides with radio emission from this region: nearly circular shell 
with a radius of 40 pc. Actually, this picture resembles that of the outer 
rims of young supernova remnants. One should realize that the life time of 
superbubbles exceeds that of supernova remnants by a large factor, result- 
ing in a particle acceleration time that is substantially longer than that of 
supernova remnants. 

Also the OB association LH9 located in the HII complex Nil in the Large 
Magel l anic C ould shows evidence for non-thermal X-ray emission |Maddox] 



et al. (2009). Interestingly, the kinetic and thermal energy of LH9 is only 



half of the energy mechanical energy transferred into the OB association 
by its stars. This is consistent with current models of CR acceleration by 
superbubbles (e.g. Bykov 20011. Overall, we would expect a galaxy with a 
high star formation rate to have a higher cosmic-ray energy density than 
a galaxy with an average star formation rate. Indeed, the starburst galaxy 

" |2009a[ ). 



NGC 253 has been detected by the HESS telescope (Acero et al. 



Also, Ackermann et al. (2011) detected the star- forming region Cygnus X in 
GeV 7-rays. 

The CR Isotope Spectrometer aboard the Advanced Composition Ex- 
plorer spacecraft found that the isotopic abundances of the CR spectrum is 
consistent with solar, with a 20% contribution of Wolf-Rayet material. Since 
Wolf-Rayet stars are stages in the evolution of OB stars, that largely reside 
in B associations and superbubbles, this supports the idea that a substantial 
fraction of the Galactic CRs are accelerated in superbubbles. 

The efficiency with which kinetic energy in superbubbles can be trans- 
ferred into CR energy are estimated to be around 20% ( Bykov |2001 ) . Super- 
bubbles are surrounded by massive shells and therefore expand much slower 
than supernova remnants reducing the adiabatic deceleration and alleviating 
the CR escape from superbubbles. Magnetic fields inside the superbubble 
should be amplified up to about 30 /iG to provide the maximal energies 
of CR protons to be about 10 s GeV and higher for heavy CR nuclei that 
will dominate the accelerated CR composition at the highest energies in the 
model. We illustrate in Figure [l6| the CR composition simulated in a model 
where CRs of energies above the knee are accelerated in a superbubble (see 
|Bykov and Toptygin[ |2001[ ). 



7.2 Reaching and passing the knee 



A decade ago, the magnetic field strength in supernova remnants was be- 
lieved to be not more than four times Galactic magnetic field as this is the 
compression factor of a non-accelerating shock. This implies a rather large 
diffusion coefficient for particles that are being accelerated and therewith a 
relative large time b efore the particles return over the shock front. [Lagagc 
and Cesarsky ( 1983 ) were the first to calculate that with these high diffusion 



coefficients, supernova remnants are not capable of accelerating particles up 
to the knee in the CR spectrum (see also section 2.3). 
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Fig. 16 CR composition abo ve the knee simulated in a m odel of CR acceleration in 
a superbubble (adapted from |Bykov and Toptygin 2001[ ). The energy dependence 
of the mean logarithm of atomic mass InA ot accelerated particles is shown. 



Bell and Lucek (2001) calculated that CRs, as they stream up and down 



the shock front, can amplify the magnetic fields themselv es by even an order 
of magnitude. This has observationally been confirmed (Vink and Laming 



[2003] ) 

nants (IVolk et al. 2005 



and has now been found for a substantial amount of supernova rem- 
and Table [3]) . This shows that from a theoretical 



perspective, supernova remnants should be able to accelerate particles up to 
the knee. 

Spectra of X-ray synchrotron emitting remnants have shown that the cut- 



off of the electron spectrum is at TeV energies (Parizot et al. 2006). This 



is still a few orders of magnitude short of the knee, but we know that the 
maximum energies of electrons are limited by the energy they lose as they 
emit synchrotron emission. Protons do not suffer from these losses and might 
well reach higher energies. 

Pion decay 7-rays are produced by protons in a more or less direct fashion. 
Therefore, observing supernova remnants with current GeV and TeV 7-ray 
telescopes would yield a wealth of information on their CR proton content. 
However, for some remnants, it is debated whether to interpret the TeV 7- 
ray data in terms of Inverse Compton scattering (which would be electrons 
again) or pion decay. In any case, even if pion decay would be the correct 
model, none of the remnants observed thus far show evidence for protons 
with energies up to 3 x 10 15 eV. 

There is, however indirect evidence for protons with these high energies in 
supernova remnants, as shown by the stripes of non-thermal X-ray emission 
on the surface of Tycho Eriksen et al. (2011 ); Bykov et al. (2011a). Moreover, 



(Ptuskin et al. 2010) found that remnants of t ype lib supernovae should be 
able to accelerate Fe ions up to 5 x 10 18 eV and Morlino and Caprioli (2011 ) 
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modeled the broadband emission of Tycho's supernova remnant and found 
that protons were at least accelerated up to 5 x 10 14 eV. 



8 Summary 

We reviewed the current observational research on the sources of Galactic 
cosmic rays. We described current shock acceleration theory, and how non- 
linear shock acceleration modifies the shock structure as well as the emission 
mechanisms that reveal the presence of cosmic rays. We use this knowledge 
to review the current observational status of cosmic-ray acceleration in su- 
pernova remnants, the Galactic center region and superbubbles. 

In these concluding remarks, we briefly point out some of the remaining 
open issues related to high Mach number shocks of supernova remnants being 
the main sources for Galactic cosmic rays. One of the big uncertainties in 
present-day cosmic-ray acceleration theory is the physics that determines 
the amount of particles getting injected into the acceleration process. This 
is particularly important, as it is one of the key ingredients for determining 
the efficiency of shocks accelerating cosmic rays. 

Also, the maximal energy to which cosmic-rays can be accelerated by 
supernova remnants is unclear. The discovery of magnetic field amplification 
makes acceleration up to, or beyond the knee, now much more plausible than 



was thought three decades ago |Lagage and Cesarsky| ( |1983| , but observation 
proof for the presence of particles with energies in excess of 10 15 eV in SNRs is 
still lacking. Third, we would like to point out that the moment at which the 
bulk of the cosmic rays escape the remnant still unclear is. In this context, 
it is useful to point out that most of the remnants detected in 7-rays are 
mature remnants, interacting with a molecular cloud, suggesting that these 
remnants still have a substantial cosmic-ray content at older ages. However, 
the maximum energies of the primary cosmic-ray particles seem in these 
mature SNRs limited to at best 1-10 GeV. 

In the near future more results are expected from the current generation 
of X-ray, and 7-ray telescopes. But future telescopes such as the Cherenkov 
Telescope Array, will greatly improve the sensitivity which gamma-rays can 
be detected. In addition, our knowledge of high Mach number shocks in 
general, and the role of cosmic rays for their thermodynamics, will improve 
with future X-ray and optical spectroscopic observations, such as with the 
Japanese X-ray telescope Astro-H and 30m class optical telescopes such as 
ELT. 

These open issues will most efficiently be solved if theory and observa- 
tions work hand-in-hand. Workshops as this in Bern, provide indispensable 
opportunities for close interaction between theorist and observers. 
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